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saturable  absorber  and  SHB  is  developed  and  studied  in  detail  both  analytically  and  numerically.  A  wide  variety  of  experimental  data  on 
multimode  regimes  is  presented.  Lasers  with  narrow  active  region  and/or  with  metal  coating  on  the  sides  tend  to  develop  a  splitting  in  the 
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1.  Introduction 


This  report  describes  in  detail  the  activity  of  Prof.  Capasso  (Harvard  University, 
responsible  for  the  experimental  part)  and  Prof.  F.X.  Kaertner  (MIT,  responsible  for  the 
modeling  part)  and  covers  the  entire  work  performed  throughout  the  duration  of  the 
U.S.  Army  Research  Office  Grant  No.  W91  INF-04-1-0253. 

Many  physical  mechanisms  can  drive  a  laser  from  a  single-mode  to  a  multi-mode 
regime.  Common  examples  are  spatial  and  spectral  hole  burning  (SHB),  saturable 
absorption,  and  self-phase  modulation  [1,2,3],  Understanding  these  mechanisms  is  of 
key  importance  to  laser  science  and  technology,  whether  one  is  interested  in 
single-mode  behavior  or  in  a  particular  multi-mode  operation  such  as  mode  locking. 
The  multimode  regimes  listed  above  are  well  understood  and  documented  both 
theoretically  and  experimentally,  in  many  types  of  lasers.  However  the  understanding 
of  multimode  regimes  in  quantum  cascade  lasers  (QCLs)  is  still  in  its  infancy,  as  these 
lasers  were  only  demonstrated  in  1994  [4],  and  studies  of  their  multimode  regimes 


commenced  even  more  recently  [5, 6, 7, 8], 

As  it  was  shown  recently  [9]  that  multimode  dynamics  in  QCLs  are  different  from  that 
of  more  common  lasers.  This  is  mainly  due  to  the  unusually  fast  gain  recovery  of  QCLs, 
which  occurs  on  a  picosecond  scale.  While  a  saturable  absorber  triggers  mode  locking 
in  lasers  with  slow  gain  recovery  (relative  to  the  roundtrip  time),  in  lasers  with  fast  gain 
recovery  a  saturable  absorber  triggers  a  mechanism  similar  to  the 
Risken-Nummedal-Graham-Haken  (RNGH)  instability  [10,11],  While  in  standard 
semiconductor  lasers  carrier  diffusion  eliminates  spatial  hole  burning,  in  QCLs  the  gain 
recovery  process  is  faster  than  carrier  diffusion,  and  spatial  hole  burning  is  dominant. 
This  paper  presents  a  detailed  study  of  multimode  regimes  in  QCLs.  The  results  of 
Ref.  [9]  are  substantiated  and  extended.  The  first  part  of  the  paper  is  theoretical  and  the 
second  one  is  experimental.  In  the  theoretical  section,  the  laser  Maxwell-Bloch 
equations  are  introduced  and  analyzed  for  both  ring  laser  cavity  and  a  Fabry-Perot 
cavity.  This  way  one  can  study  the  interplay  of  coherent  phenomena  and  spatial  hole 
burning.  A  saturable  absorber  is  added  to  the  model  as  well.  The  model  is  studied 
analytically  and  numerically.  The  stability  region  of  a  continuous  wave  (CW)  solution 
is  found.  It  is  shown  that  in  a  ring  cavity,  the  presence  of  a  saturable  absorber  lowers  the 
threshold  of  the  RNGH  instability  from  about  nine  times  above  laser  threshold  to 
arbitrarily  low  above  laser  threshold,  depending  on  the  strength  of  the  absorber. 
However,  the  nature  of  the  instability  remains  the  same:  the  population  inversion 
begins  to  oscillate  at  the  Rabi  frequency,  modulating  the  gain  in  the  laser.  The  result  is 
sidebands  around  the  original  CW  mode,  separated  from  it  by  roughly  the  Rabi 
frequency.  A  Rabi  splitting  in  the  spectrum  is  the  primary  signature  of  the  RNGH 
instability. 

It  is  then  shown  that  the  ring-cavity  Maxwell-Bloch  model  with  a  saturable  absorber 
can  explain  the  Rabi  splitting  and  the  lowering  of  the  threshold,  but  cannot  explain  a 
key  feature  in  the  experimental  spectra.  The  latter  exhibit  two  relatively  equal  groups  of 
modes  with  a  gap  in  between,  whereas  the  ring-cavity  Maxwell-Bloch  model  predicts  a 
large  central  mode  in  the  spectrum  with  two  sidebands.  In  view  of  this  discrepancy,  the 
Maxwell-Bloch  equations  were  extended  to  include  coupling  between 
counter-propagating  modes  in  a  Fabry-Perot  cavity,  which  is  of  course  a  more  adequate 
model  for  QCLs.  This  model  takes  into  account  the  development  of  SHB.  The  inclusion 
of  SHB  generates  theoretical  spectra  with  a  Rabi  splitting  and  without  a  central  peak  in 
the  spectrum,  in  agreement  with  experiments. 

After  presenting  the  theoretical  results,  a  thorough  study  of  the  experimental 
phenomenology  is  presented.  It  is  shown  that  QCLs  with  narrower  active  regions  tend 
to  exhibit  a  more  pronounced  Rabi  splitting  than  lasers  with  wider  active  regions. 
Lasers  with  a  wider  active  region  tend  to  exhibit  multimode  spectra  that  are  governed 
by  SHB.  The  explanation  we  propose  to  this  behavior  is  that  narrow  QCLs  have  a 
stronger  saturable  absorption  effect  than  wider  ones  as  a  result  of  a  Kerr-lensing  effect 
[3,5],  In  this  case  a  nonlinear  index  enhancement  in  the  waveguide  core  gives  rise  to  an 
increased  overlap  of  the  transverse  laser  mode  with  an  active  region  and  a  reduced 
overlap  with  lossy  sidewalls,  leading  to  an  additional  enhancement  of  the  saturable 
absorption.  The  idea  of  saturable  absorption  by  Kerr-lensing  is  further  supported  by  the 
fact  that  QCLs  with  metal  coating  on  the  sides  of  the  ridge  have  a  stronger  RNGH 
behavior  than  lasers  without  metal  coating:  The  metal  coating  enhances  saturable  losses 
originated  from  Kerr-lensing. 

In  the  last  part  of  the  paper,  the  temperature  dependence  of  the  multimode  behavior  in 
QCLs  is  studied  experimentally.  It  is  found  that  at  higher  temperature  the  multimode 
behavior  tends  to  be  governed  by  the  RNGH  instability,  whereas  at  lower  temperatures 


it  is  governed  by  SHB.  This  behavior  can  be  explained  by  the  fact  that  at  higher 
temperatures  hot  carriers  populate  more  states  in  the  injector  superlattice,  creating 
aditional  quasi-resonant  absorption  transitions  between  ground  and  excited  minibands. 
This  leads  to  an  additional  saturable  absorption  of  laser  radiation. 

This  paper  is  organized  as  follows:  Section  2  gives  a  brief  survey  of  prior  work  on 
multimode  regimes  in  QCLs,  Sections  3-5  study  theoretically  the  Maxwell-Bloch 
equations  in  a  ring  cavity  and  a  Fabry-Perot  cavity,  Sections  6-8  summarize  the 
experimental  study,  and  Section  9  is  a  brief  summary. 

Note  that  most  of  the  Quantum  Cascade  Lasers  (QCLs)  investigated  during  this  project 
were  grown  by  Dave  Bour  of  Agilent  Technologies.  This  fruitful  collaboration  with 
industry  has  led  to  the  publication  of  several  articles  (see  for  example  References  [28] 
and  [29]  and  the  papers  attached  after  Appendix  A  and  B).  The  wafers  were  grown  by 
metal  organic  vapor  phase  epitaxy  (MOVPE),  a  crystal  growth  technique  fitting  the 
requirements  of  large-scale  production.  The  material  quality  and  the  performance  in 
continuous  mode  were  shown  to  be  comparable  to  or  better  than  that  of  state-of-the-art 
QCLs  grown  by  the  standard  deposition  technique  for  these  devices,  namely  molecular 
beam  epitaxy  or  MBE. 


2.  Prior  work  on  multimode  regimes  in  QCLs 

Multimode  regimes  in  QCLs  were  observed  in  a  series  of  recent  works  [5, 6, 7, 8],  In 
Ref.  [5],  for  example,  it  was  observed  that  at  a  certain  pumping  current  above  lasing 
threshold,  QCLs  cease  to  operate  in  CW  and  develop  a  multimode  regime.  This 
multimode  regime  was  characterized  by  a  broadband  optical  spectrum  and  a  narrow 
(less  than  100  kHz)  radio  frequency  (RF)  beat  note  in  the  power  spectrum. 

The  narrow  beatnote,  whose  width  is  1/105  of  the  central  frequency,  shows  that  the 
waveform  of  the  electric  field  circulating  in  the  laser  cavity  was  stable  over 
approximately  105  roundtrips.  In  other  words,  the  phase  relationships  between  the 
longitudinal  modes  were  stable  for  about  105  roundtrips.  The  modes  were  therefore 
locked.  However  in  order  to  characterize  the  waveform  which  is  circulating  in  the  laser 
and  to  see  if  it  is  indeed  an  isolated  pulse,  as  in  traditional  mode  locking,  one  has  to 
apply  pulse  characterization  techniques,  such  as  second  order  autocorrelation. 

At  the  time  when  the  experiment  in  Ref.  [5]  was  performed,  no  second-order 
autocorrelation  apparatus  was  available.  However  second  harmonic  generation  from 
QCLs  [6]  provide  some  information  for  pulse  characterization.  This  measurement 
shows  an  increase  by  more  than  a  factor  of  5  in  the  second  harmonic  signal  as  the 
multimode  behavior  sets  in.  This  increase  indicates  that  the  duty  cycle  of  the  pulses  was 
roughly  5.  Since  the  number  of  modes  in  the  spectrum  is  more  than  5,  we  conclude  that 
not  all  the  modes  were  locked  into  a  pulse.  Without  better  pulse  characterization  data, 
one  can  infer  that  the  laser  could  have  had  more  than  one  pulse  per  roundtrip. 

It  should  also  be  noted  that  traditional  mode  locking,  with  a  single  pulse  per  roundtrip, 
requires  that  the  gain  recovery  time  be  longer  than  the  cavity  roundtrip  [3].  In  QCLs 
this  condition  is  violated,  and  therefore  according  to  the  existing  theory  one  cannot 
expect  mode  locking  with  a  single  pulse  per  roundtrip  in  QCLs.  One  can  expect 
multiple  pulses  per  roundtrip. 

In  view  of  the  above,  one  can  see  that  the  nature  of  various  multimode  regimes  in  QCLs 
requires  further  elucidation.  This  is  the  main  goal  of  the  present  work. 


3.  Maxwell-Bloch  equations  in  a  Fabry-Perot  cavity 


In  this  section  we  derive  the  Maxwell-Bloch  equations  in  a  Fabry-Perot  cavity.  We 
model  the  gain  medium  of  QCLs  as  a  two  level  system,  described  by  the  Bloch 
equations 
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where  pab  is  the  off-diagonal  element  of  the  density  matrix,  A=pbb-paa  is  the  population 
inversion,  co  and  p  are  the  resonant  frequency  and  the  dipole  matrix  element  of  the 
lasing  transition,  Ti  and  T2  are  the  longitudinal  and  transverse  relaxation  times,  and  Ap 
is  the  steady-state  inversion  at  E=0  which  characterizes  the  pumping  rate.  The  last  term 
in  Eq.  (2)  is  added  phenomenologically  and  accounts  for  spatial  diffusion  of  the 
inversion  due  to  carrier  diffusion.  D  is  the  diffusion  coefficient.  E  is  the  electric  field, 
which  is  assumed  to  satisfy  the  wave  equation 
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N  is  the  number  of  two-level  systems  per  unit  volume,  which  equals  to  the  doping 
density  multiplies  the  number  of  QCL  stages.  T  is  the  overlap  factor  between  the  optical 
mode  and  the  active  region,  and  n  is  the  background  refractive  index. 

Eqs.  (1-2)  describe  an  open  two-level  system  [41]:  The  total  number  of  electrons  in  the 
system  is  not  conserved,  but  rather  can  flow  in  and  out  of  the  system.  In  general  the  two 
levels,  the  upper  one  b  and  the  lower  one  a,  can  have  different  relaxation  times.  In  this 
situation  Eq.  (2)  should  be  replaced  by  two  equations,  each  with  its  own  relaxation 
times.  In  QCLs  indeed  the  upper  and  lower  levels  have  different  relaxation  times. 
However  we  neglect  this  difference  for  the  sake  of  simplicity  and  assign  to  them  the 
same  value  Ti.  The  generalization  of  the  model  (1-2)  to  a  model  with  two  different 
relaxation  time  is  straightforward,  as  also  is  the  generalization  to  a  model  with  more 
than  two  levels. 

We  now  make  the  following  set  of  Ansatz: 


E(z,t)  =  ^ 

1 


+E+(z.t)e' 


+ 

i  i) 


Pab(z,t)  =  j]+(zJ)eKu>t-k!:) +  i]-(z,t)ei^t+kz)  (5) 
A(*,i)  =  Ao(M)  +  Aa(«,0e!M**  +  (6) 


where  k  =  co  n  /c.  E±,  r\±,  Ao  and  A2  are  assumed  to  vary  slowly  varying  in  time  and  space, 
on  the  scale  defined  by  co  and  k  respectively.  The  quantities  with  a  +  (-)  subscript 
represent  waves  traveling  to  the  positive  (negative)  z  direction.  Eq.  (7)  allows  taking 
SHB  into  account,  A2  being  the  envelope  of  the  inversion  grating.  Note  that  Eq.  (7)  can 
be  extended  by  adding  terms  proportional  to  e4,kz ,  e6,h  etc.  Neglecting  these  terms 
means  that  higher  spatial  frequencies  on  the  inversion  are  neglected.  Due  to  the 
nonlinearity  of  the  gain  medium  such  frequencies  can  appear  if  the  gain  is  heavily 
saturated,  but  they  are  neglected  in  our  analysis  for  the  sake  of  simplicity  and  since  the 
pumping  in  our  system  is  never  much  higher  than  the  laser  threshold. 


Substituting  Eqs.  (4-7)  into  Eqs.  (1-3)  and  making  the  slowly- varying  envelope 

approximation,  we  obtain  the  following  set  of  equations: 

n .  kNuV  1  _  _  .  _ 

-OtE±  =  ^dzE±-  i- - Tt}±  -  -i{E+,E_)E±  (7) 

c  2 

=  ^(Ao£±  +  A 2=  £t)  -  (8) 

2n  i  2 

a,i0  =  '  A"  +  L(£>+  +  £!>)-  -  c.=.)  (9) 

a, At  =  ±i|(£;>, -  4i-2CAj(io) 

where  we  have  introduced  the  notation  A+2  =  A2 ,  A“  =  A*,  in  order  that  Eqs.  (8-11)  can 
be  written  more  compactly.  The  last  term  in  Eq.  (8)  has  been  added,  and  represents  loss. 
The  loss  [  is  allowed  to  depend  on  the  field  to  represent  phenomena  such  as  optical 
saturation. 

In  QCLs  the  laser  cavity  is  formed  by  the  two  interfaces  between  the  semiconductor 
medium  and  the  outside  air,  one  located  at  z=0,  and  the  other  one  at  z=L.  L  is  the  length 
of  the  cavity.  At  each  interface  the  Fresnel  reflection  law  dictates  the  following 
relations: 


E_(LJ] 

n  —  1 

)=  £+(£■*) 
n  +  1 

(11) 

E+(0.t 

n  —  1 

)  = - — £_(0.t) 

+  1 

(12) 

In  what  follows  we  study  analytically  and  numerically  the  model  introduced  in  this 
section.  We  begin  with  a  simplified  case,  namely  the  standard  Maxwell-Bloch  equation 
in  a  ring  cavity.  In  addition  to  briefly  reviewing  known  results  about  the  RNGH 
instability,  we  study  the  effect  of  a  saturable  absorber  on  the  latter. 

4.  Ring  cavity 

In  this  section  we  consider  a  ring  cavity,  where  SHB  does  not  exist  because  standing 
waves  cannot  form.  The  aim  is  to  understand  the  interplay  of  coherent  effects  and  a 
saturable  absorber  alone,  before  further  complications  due  to  SHB  come  in.  We  shall 
see  that  without  SHB,  the  qualitative  agreement  between  theory  and  the  experiments  on 
QCLs  is  not  complete.  After  introducing  SHB  in  the  next  section,  the  agreement  is 
much  more  satisfactory. 

4.1.  The  RNGH  instability  with  a  saturable  absorber 

Dropping  all  the  quantities  with  a  subscript  from  Eqs.  (8-11),  one  arrives  at  the 
standard  Maxwell-Bloch  equations,  with  a  saturable  absorber  added. 

-0,E  =  -d,E  -  -  t(<0  -  7|E|2)E  (13) 

dm  =  (14) 

2  n  1 2 

dt A  -  ^  +  t( E*ij-c.c.).  (15) 

U 

The  saturable  absorber  is  approximated  to  the  lowest  order  in  E,  and  is  characterized  by 
y,  which  is  often  referred  to  as  the  self  amplitude  modulation  coefficient,  /o  is  the  linear 


loss.  Ath  is  the  lasing  threshold  value  of  Ap  for  y  =0,  given  by: 


A-Ay.-fT? 


( 16) 


g(Q)  for  p— 2  and  the  parameters  in  Table  I,  apart  fromy.  The  latter  is  y— 0  (solid), 
y  =  10 ~9cm/V2  (dashed)  and  y  =  2xl0 ~9cm/V2  (dotted). 


Linear  stability  analysis  of  Eqs.  (14-16)  can  be  carried  out.  The  gain  of  a  perturbation  at 
the  frequency  Q  (relative  to  the  resonance  frequency  co)  is  approximately  given  by 

(QTi  +  i)f>r2  -  2(p  -  1) 


g(Q)  =  -  —  Re 


+ 


(f!Xi  +  i)(£2T2  +i)  —  (p  —  1) 
'fH2(p  -  1)  (flTj  +  i)(3QX2  +  2i)  -  4(p  -  1) 
pi^'T 1T2  (t2Xi  +  i)  (£TX2  +  i)  —  p  +  1 


(17) 


The  approximations  made  in  the  derivation  of  Eq.  (18)  mainly  include  assuming  that 
the  photon  lifetime  in  the  empty  cavity  is  much  longer  than  Ti  and  T2.  This 
approximation  is  excellent  for  QCLs.  p  is  the  pumping  above  lasing  threshold  (for  y=0): 


P  = 


Ath* 


(18) 


Figure  1  shows  g(Q)  for  different  values  of  y  at  p=2.  The  parameters  used  in  all 
examples  in  this  paper,  unless  specified  otherwise,  are  given  in  Table  I.  The  effect  of  y 
is  to  increase  g(Q)  more  or  less  uniformly  across  the  frequency  domain.  In  particular,  it 
can  bring  g(Q)  above  zero,  thereby  triggering  an  instability,  even  when  the  laser  is  not 
pumped  as  high  above  threshold.  The  reason  why  a  saturable  absorber  lowers  the 
RNGH  threshold  is  that  a  saturable  absorber  itself  always  favors  a  multimode  regime  to 
a  single  mode  one.  It  introduces  a  frequency-independent  parametric  gain.  The  latter  is 
added  to  the  RNGH  parametric  gain  from  Fig.  1,  bringing  it  above  threshold.  Note, 
however,  that  the  instability  still  starts  by  developing  Rabi  sidebands  around  the  CW 
lasing  frequency.  In  this  sense,  it  can  be  interpreted  as  a  modified  version  of  the  original 
RNGH  instability,  rather  than  as  a  modulation  instability  caused  by  the  saturable 
absorber  alone. 


Quantity 

Symbol 

Value 

Gain  recovery  time 

7i 

0.5  ps 

Dephasing  time 

t2 

0.0+57  ps 

Linear  cavity  loss 

5  cm-1 

Transition  dipule  element 

jU 

2.54  rnnx  e 

Background  refractive  index 

n 

3 

Cavity  length 

L 

3  mm 

Saturable  absorber  coefficient 

7 

1.  The  parameters  used  in  all  calculations  and  simulations  in  this  paper,  unless  indicated 

otherwise. 

The  dependence  of  the  threshold  for  instability  on  y  is  shown  in  Fig.  2.  For  y=0  we 
recover  the  standard  RNGH  instability,  which  occurs  at  the  pumping  of  slightly  above  9 
times  above  lasing  threshold.  For  y— *  oo  the  instability  threshold  approaches  the  lasing 
threshold. 


2.  The  pumping  ratio  p  at  which  the  RNGH  instability  sets  in  as  function  of  the  saturable 
absorber  coefficient  y.  Solid  line  -  exact  result  from  linear  stability  analysis.  Dashed  line  - 
approximate  result  (Eq.  (23)).  The  parameters,  apart  from  y,  are  given  in  Table  I. 

The  results  in  Fig.  1  and  Fig.  2  were  obtained  from  numerical  solutions  of  the  algebraic 
equations  involved  in  the  stability  analysis.  Since  even  the  approximate  expression  of 
g(Q)  (Eq.  (18))  is  not  very  simple,  it  is  useful  to  derive  some  approximate  simple 
expressions  for  the  properties  of  the  instability.  Here  we  only  give  the  results. 


3.  Qmax  as  function  of  the  pumping.  Thick  line:  exact,  dotted  line:  Eq.  (21).  The  parameters  are 
from  Table  I  with  the  only  difference  that  y=0 


g(Q)  (e.  g.  in  Fig.  1)  has  a  local  minimum  at  Q=0,  and  peaks  at  approximately 


Lmaj: 


(19) 


=  0 


' 'Rabi 


J 

V  p- 1 


The  position  of  the  peak  depends  weakly  on  y  within  the  parameter  range  of  interest  to 
our  system.  Note  that  for  y=0  one  has 


P-1 


ry  _  H _ 

Rabi  -  y  jj'  ■ 


(20) 


The  dependence  of  |Qmax|on  p  is  shown  in  Fig.  3. 

The  instability  threshold  is  approximately  given  by 

(21) 

Figure  2  shows  that  Eq.  (23)  fairly  well  approximates  the  exact  threshold  condition. 


4.2.  Numerical  results 

The  linear  stability  analysis  given  above  can  only  find  the  threshold  condition,  and  does 
not  describe  the  behavior  of  the  laser  above  the  point  of  instability.  Here  we  therefore 
present  the  results  of  a  numerical  solution  of  Eqs.  (14-16). 

According  to  Ref.  [30],  the  state  of  the  laser  can  change  either  discontinuously  or 
continuously  at  the  RNGH  instability  threshold.  In  the  language  of  phase  transitions, 
this  would  be  a  first  or  second  order  phase  transition  respectively.  The  criterion  that 
determines  the  order  of  the  transition  involves  the  laser  cavity  length,  but  has  not  been 
given  explicitly  in  the  literature  known  to  us. 

According  to  the  numerical  simulations,  for  the  parameters  corresponding  to  our  QCLs, 
including  the  saturable  absorber,  the  transition  is  continuous  (second  order).  Below  the 
instability  threshold  the  lasing  is  CW,  and  as  the  threshold  is  crossed,  the  Rabi 
sidebands  around  the  central  CW  mode  grow  continuously.  This  behavior  is 
demonstrated  in  Fig.4.  Note  that  standard  mode  locking,  with  a  saturable  absorber  and 
slow  gain  recovery,  is  a  discontinuous  (first-order)  transition  [43,44], 

The  spectra  in  Fig.4  have  two  groups  of  modes  separated  by  roughly  twice  the  Rabi 
frequency.  In  this  respect  they  resemble  the  experimental  spectra.  However  they  have  a 
strong  CW  peak  in  between,  a  feature  which  is  not  shared  by  the  experimental  spectra. 
As  we  show  in  the  next  section,  in  a  Fabry-Perot  cavity  with  SHB,  the  central  CW  peak 
disappears. 


4.  Spectra  obtained  from  a  numerical  solution  of  the  Maxwell-Bloch  Eqs.  (13-15).  The  parameters 

are  given  in  Table  I. 


5.  Fabry-Perot  cavity 

Spatial  hole  burning  is  associated  with  A2  in  Eq.  (11).  Intuitively,  A2  is  the  amplitude  of 
the  grating  that  couples  the  two  propagation  directions  in  the  laser.  The  parameter  that 
controls  the  strength  of  SHB  is  D:  in  the  limit  of  D— >oo,  A2  approaches  zero.  In  order  to 
better  understand  the  interplay  between  SHB  and  the  RNGH  instability,  we  present 
now  the  results  of  analytical  and  numerical  studies  of  Eqs.  (8-11).  We  start  with  linear 
stability  analysis  and  here  we  only  give  the  results. 

Before  proceeding,  we  define  the  lifetime  of  the  gain  grating  Tg  as  seen  in  Eq.  11,  as 

T"1  =  Tf1  +  ±k‘2D. 

Tg  is  the  parameter  that  determines  the  strength  of  spatial  hole  burning.  Tg  can  therefore 
range  from  zero  (no  SHB)  to  Ti  (strongest  of  SHB).  The  diffusion  coefficient  D  can  be 
estimated  from  the  Einstein  relation.  In  AlInAs-InGaAs  heterostructures,  the  electron 
mobility  can  reach  as  high  as  16000  cm  /sec/V  at  room  temperature  [45].  Thus  one  has 
400  cm2/sec,  k=3.7*104  cm'1,  which  roughly  corresponds  to  a  vacuum  wavelength  of  5 
pm,  we  obtain  4k2D~  2  THz.  With  7j_1  *  277/z  we  find  that  Tg~0.25  ps.  Note  that  the 
mobility  used  here  was  relatively  high,  the  wavelength  was  on  the  short  side  of  the  scale 
and  room  temperature  was  taken.  Therefore  in  reality  Tg  is  closer  to  Ti.  It  therefore 
follows  that  due  to  the  fast  gain  recovery  of  QCLs,  carrier  diffusion  does  not  eliminate 
spatial  hole  burning.  This  is  in  contrast  with  diode  lasers. 


5.1.  Linear  stability  analysis 


Linear  stability  analysis  of  Eqs.  (8-11)  gives  two  families  of  unstable  modes.  One  is 
associated  with  the  RNGH  instability,  and  the  other  one  -  with  spatial  hole  burning. 
The  first  family  is  very  similar  to  the  case  of  a  ring  cavity  studied  earlier,  with  small 
differences.  The  second  one  can  also  be  derived  and  the  result  is  shown  in  Fig.  5. 


5.  g(f2)  of  the  perturbation,  associated  with  spatial  hole  burning,  p  and  Tg  corresponding  to 

every  curve  is  denoted  in  the  legend. 


Figure  4  shows  the  gain  of  a  perturbation  around  a  CW  solution.  The  shorter  is  Tg  and 
the  smaller  is  the  pumping,  the  smaller  is  the  gain  of  the  instability.  For  parameters 
typical  to  our  QCLs,  the  SHB  instability  occurs  a  few  percents  above  lasing  threshold. 
The  gain  curve  in  Fig.  5  exhibits  two  peaks  and  a  dip  at  Q=0.  For  T2«Ti  and  p- 1 « 1 , 
the  location  of  the  peaks  is  given  by: 
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_L  / IEL 

Ti  V  37Y74 ' 


(22) 


Note  that  Eq.  (22)  gives  a  smaller  frequency  than  Eq.  (19).  In  addition,  the  splitting  in 
Fig.  5  scales  like  the  square  root  of  the  Rabi  frequency.  CW  is  destabilized  when  the 
cavity  admits  a  mode  for  which  g(Q)  in  Fig.  5  is  positive. 

5.2.  Numerical  results 


Wavenu  mber  (cm1) 


6.  Results  of  numerical  simulations  of  the  spectra  based  on  the  Maxwell-Bloch  equations  including 
a  saturable  absorber  and  spatial  hole  burning  for  different  values  of  the  current  density 

normalized  to  the  threshold  value. 

The  moment  CW  is  destabilized,  studying  Eqs.  (7-10)  requires  a  numerical  simulation. 
The  results  of  such  a  simulation  are  shown  in  Fig.  5.  The  parameters  are  given  in  Table 
I,  with  the  only  difference  that  y=0  was  used. 

Fig.  6  shows  a  clear  pattern  in  the  spectrum.  This  pattern  appears  only  after  very  long 
averaging  (order  of  a  microsecond)  of  the  spectrum.  Such  averaging  is  appropriate, 
since  similar  averaging  occurs  when  the  spectra  are  measured.  The  pattern  depends  on 
T2  and  Co,  but  very  weakly  depends  on  Ti.  A  similar  pattern  occurs  in  the  experimental 
spectra.  However  we  were  not  able  to  understand  its  origin. 

5.3.  Saturable  absorber 

The  combined  effect  of  SHB  and  a  saturable  absorber  is  demonstrated  in  Fig.  7.  One 
can  observe  two  effects.  First  is  spectral  broadening  due  to  SHB,  similarly  to  Fig.  6. 
Second  however  is  the  appearance  of  a  splitting  in  the  spectrum.  This  splitting  is 
roughly  equal  to  twice  the  Rabi  frequency.  Note  however  that  contrary  to  a  ring  laser,  in 
a  Fabry-Perot  laser  "the  Rabi  frequency"  is  not  a  well  defined  concept.  Since  a  standing 
wave  is  formed  in  the  cavity,  the  electric  field  and  thus  the  Rabi  frequency  depend  on 
the  position  in  the  cavity.  This  dependence  is  even  stronger  when  the  the  mirrors  have  a 
relatively  low  reflection  coefficient,  since  the  field  amplitude  even  more  strongly 
depends  on  position. 

In  the  previous  section  we  saw  that  a  saturable  absorber  lowers  the  threshold  of  the 
RNGH  instability.  Here  we  see  that  SHB  suppresses  the  central  peak  seen  in  Fig.4  and 


replaces  it  by  a  minimum  in  the  spectrum.  In  the  next  sections  we  see  that  Fig.  7  agrees 
well  with  experimentally  measured  spectra,  at  least  for  devices  where  the  RNGH 
behavior  was  dominant. 
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7.  Results  of  numerical  simulations  of  the  spectra  based  on  the  Maxwell-Bloch  equations  including 
a  saturable  absorber  and  spatial  hole  burning  for  different  values  of  the  current  density 

normalized  to  the  threshold  value. 

6.  Experimental  results:  ridge  lasers 

We  first  study  the  multimode  regimes  in  standard  ridge  QCLs,  in  which  the  sidewalls  of 
the  laser  ridges  are  covered  by  a  thick  layer  of  electrically  plated  gold  contact.  This  acts 
as  a  Kerr-lens  type  saturable  absorber;  see  Fig.  8  for  a  typical  laser  cross  section.  The 
Kerr-lensing  mechanism  is  explained  in  detail  in  Appendix  B.  The  active  region  of  the 
samples  tested  is  based  on  a  three-quantum-well  design  emitting  at  a  wavelength 
X~8pm  [31].  The  wafer  was  grown  by  metalorganic  vapor  phase  epitaxy  (MOVPE). 
Fig.  9  a  shows  the  voltage-current  (V-I)  and  light-current  (L-I)  characteristics  of  a  10 
pm  wide  laser  operated  in  CW  at  200K,  and  b  shows  the  corresponding  optical  spectra. 
The  laser  was  cleaved  into  a  2  mm  long  bar  and  soldered  with  Indium  onto  a  copper 
heat  sink.  The  optical  power  was  measured  by  an  OPHIR  thermal  head  powermeter 
with  a  collection  efficiency  of  nearly  100%.  The  spectra  were  measured  by  a  Nicolet 
Fourier  transform  infrared  spectrometer  (FTIR)  equipped  with  a  deuterated  triglycine 
sulphate  (DTGS)  detector. 


8.  Scanning  electron  microscope  (SEM)  image  of  the  cross  section  of  a  10  [Trial  mode]m  wide 

ridge  QC  laser. 


9.  (a)  V-I  and  LI  curves  and  (b)  optical  spectra  vs.  pumping  ratio  Q'/jth)  above 
threshold  obtained  in  CW  at  200K  of  a  10  pm  wide  ridge  laser  emitting  at  8pm.  (c) 
Spectral  splitting  and  twice  the  Rabi  frequency  vs.  square  root  of  output  power 
collected  from  a  single  facet.  The  different  quantities  reported  on  the  graph  were 
deduced  from  the  experimental  data  show  in  (b).  The  dashed  line  is  the 
least-square  linear  fit  of  the  data. 


As  shown  in  Fig.  9  b,  the  laser  spectrum  is  single  mode  close  to  laser  threshold.  It 
broadens  and  splits  into  two  separated  humps  as  the  pumping  current  increases.  The 
separation  between  the  two  peaks  of  the  two  humps  increases  linearly  with  the  square 
root  of  the  collected  output  power  from  one  facet,  as  shown  in  Fig.  9  c.  The  Rabi 
angular  frequency  can  be  calculated  from  the  collected  output  power,  using  the  formula 
Q  Rabi  =  juE/h,  where  p  is  the  electron  charge  times  the  matrix  element  of  the  laser 

transition  (=1.9  nm  for  this  particular  device).  Iave  is  the  average  intracavity  intensity  in 
the  gain  region,  which  can  be  derived  from  the  measured  output  power  [32],  For  all 
values  of  intensity  corresponding  to  the  spectra  reported  in  Fig.  9  b,  ClRabi  /  In  was 

calculated,  multiplied  by  a  factor  of  two  and  then  added  to  Fig.  9  c  (solid  line). 
Reasonably  good  agreement  is  found  between  the  experimental  splitting  and  twice  the 
estimated  Rabi  frequency.  The  error  bars  of  the  spectrum  splittings  come  from  the 
uncertainty  in  determining  the  exact  position  of  the  peaks,  which  is  the 


full-width-at-half-maximum  (FWHM)  of  the  humps.  As  mentioned  previously  in  the 
theoretical  section,  the  RNGH  instability  predicts  that  large  intracavity  intensity  will 
result  in  parametric  gain  at  frequencies  detuned  from  the  maximum  of  the  gain  curve  by 
the  Rabi  frequency.  The  measured  spectra  thus  show  strong  indication  of  the  RNGH 
instability  in  ridge  QCLs. 


9.  (a)  V-I  and  LI  curves  and  (b)  optical  spectra  vs.  pumping  ratio  (j/jth)  above 
obtained  from  a  10  pm  wide  ridge  QCL  operated  in  CW  at  77K  emitting  at  8  pm 
Note  that  mode  hoping  occurs  while  there  is  a  jump  in  the  LI  curve. 

The  lowering  of  the  RNGH  instability  threshold  in  our  QCLs  is  due  to  the  presence  of  a 
saturable  absorber.  This  phenomenon  is  demonstrated  analytically  in  the  theoretical 
section.  Such  a  saturable  absorption  mechanism  in  our  experiments  is  likely  to  come 
from  Kerr-lensing,  caused  by  a  nonlinear  (i.e.  intensity  dependent)  refractive  index  n2l 
in  the  active  region  [5],  As  the  light  intensity  increases,  the  mode  becomes  more 
confined  in  the  plane  transverse  to  the  propagation  direction,  and  the  net  gain  it 
undergoes  also  increases.  The  reason  is  twofold:  First,  the  mode  overlaps  more  with  the 
active  region,  leading  to  a  larger  modal  gain  (this  mechanism  is  often  called  "soft 
Kerr-lensing").  Second,  the  overlap  with  the  metal  contacts  is  reduced,  leading  to 
smaller  losses.  A  detailed  analysis  is  given  in  Appendix  B. 


11.  The  experimental  setup  of  a  two-photon  autocorrelation  measurement.  (Inset)  Conduction 
band  diagram  of  the  two-photon  QWIP  showing  three  equidistant  energy  levels. 

The  same  RNGH  splitting  in  spectra  is  observed  in  many  different  devices,  from  wafers 
grown  by  both  molecular  beam  epitaxy  (MBE)  and  MOVPE.  Fig.  10  shows  the  spectra 


of  a  laser  fabricated  from  a  MBE-grown  wafer  with  the  same  active  region  design  as  in 
Fig.  7,  taken  in  continuous  wave  at  77K.  The  spectrum  starts  from  single  mode  close  to 
the  laser  threshold,  and  the  mode  hops  at  a  pumping  ratio  of  j/jth  =1 .66,  where  there  is  a 
correponding  kink  in  the  L-I  curve  at  330mA.  At  590mA,  there  is  another  kink  in  the 
L-I  curve,  relating  to  the  broadening  of  the  spectrum  and  a  hop  in  the  mode  center. 
Mode-hopping  and  its  relationship  to  kinks  in  the  L-I  curves  are  common  phenomena  in 
semiconductor  lasers  [33,34,35],  but  its  cause  in  QC  lasers  is  not  yet  studied  and  is 
beyond  the  focus  of  this  paper.  After  the  spectrum  broadens,  it  forms  two  separated 
humps  whose  peaks  shift  apart  with  increased  optical  power,  similar  to  the 
MOVPE-grown  case. 

Characterization  of  our  ultra-short  pulses  was  done  using  the  standard  method  of 
second-harmonic  interferometric  autocorrelation.  The  setup  is  based  on  a  Michelson 
interferometer  in  which  the  input  beam  is  split  into  two  and  one  of  them  is  delayed  by 
x.Once  recombined,  the  two  pulses  are  sent  colinearly  first  into  a  nonlinear  crystal,  and 
then  a  fdter,  which  allows  only  the  second-harmonic  generation  (SHG)  component  to 
be  detected.  One  can  test  if  it  is  an  isolated  pulse  from  the  ratio  between  the  interference 
maximum  and  the  background  (see  Appendix  A).  The  pulse  duration  can  also  be 
determined.  However,  due  to  the  extremely  low  SHG  conversion  efficiency  of  mid-IR 
in  nonlinear  crystals,  the  conventional  setup  is  not  feasible.  To  overcome  this  problem, 
we  use  a  two-photon  quantum  well  infrared  photodetector  (QWIP)  which  converts  the 
second-harmonic  signal  electrically  [36,37]  instead  of  using  a  nonlinear  crystal  plus  a 
linear  detector.  The  energy  diagram  of  one  period  of  the  multi  quantum  well  detector 
under  bias  is  shown  in  the  inset  of  Fig.  11.  The  first  three  electronic  states  are  nearly 
equidistant  in  energy.  When  electrons  in  the  doped  quantum  wells  absorb  two  photons 
simultaneously  and  the  detector  is  biased  (1-3  V),  a  photocurrent  is  generated  and  the 
signal  can  be  detected  by  use  of  a  pre-amplifier  and  a  lock-in  amplifier.  This 
experimental  setup  is  diagramed  in  Fig.  11.  The  second-order  autocorrelation  trace  of 
the  MBE-grown  device  mentioned  above  is  shown  in  Fig.  12.  Interference  fringes  are 
observed  when  the  delay  time  from  one  arm  of  the  autocorrelator  is  equal  to  the 
multiples  of  the  cavity  roundtrip  time.  The  ratio  between  the  maximum  of  the 
interference  fringes  and  the  background  is  smaller  than  8,  and  the  autocorrelation  trace 
has  some  features  between  the  cavity  roundtrip  times,  indicating  that  the  multimode 
regime  observed  in  this  device  is  not  stable  mode  locking  with  a  single  pulse  per 
roundtrip.  The  meaning  of  the  autocorrelation  traces  are  discussed  in  Appendix  A,  and 
simulated  traces  are  given  therein. 


12.  A  second-order  autocorrelation  trace  of  a  8  flm  wavelength  ridge  QC  laser  (wafer  #  2743) 
under  the  condition  of  RNGH  instability.  (Inset)  Microwave  spectrum  of  photocurrent  generated 
by  a  similar  laser  (wafer  #  2721)  under  the  condition  of  RNGH  instability  (measured  with  a  68  kHz 

resolution  bandwidth). 


In  addition  to  second-order  autocorrelation,  the  microwave  spectrum  of  the  laser  output 
was  also  measured  with  an  ultrafast  QWIP  [38]  whose  bandwidth  is  52  GHz.  The  laser 
output  is  sent  directly  to  the  ultrafast  QWIP,  and  the  resulting  photocurrentis  displayed 
in  a  spectrum  analyzer.  Fig.  12  shows  the  beat  note  signal  of  the  MOVPE-grown 
sample  at  pump  current  800  mA  at  77K.  A  steady  peak  with  FWHM  of  13  MHz  at 
22.01  GHz,  which  corresponds  to  the  cavity  round  trip  frequency  of  the  2  mm  long 
laser  (background  refractive  index  n=3),  is  observed  on  the  spectrum  analyzer.  It 
indicates  a  modulation  of  the  laser  output  at  the  cavity  roundtrip  frequency,  and  thus  at 
least  partial  phase  locking  between  the  longitudinal  modes:  The  phase  relationships 
between  the  modes  are  stable  for  about  103  roundtrips. 
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13.  Optical  spectra  vs.  pumping  ratio  ([Trial  mode])  above  threshold  obtained  in  CW  at  77K 
with  a  15  pm  wide  ridge  laser  emitting  at  8  pm  (wafer#  2721). 

As  the  ridge  width  is  increased,  the  RNGH  instability  and  the  Rabi  splitting  in  the 
spectrum  are  suppressed.  Fig.  13  shows  the  spectra  of  a  device  processed  from  the  same 
wafer  as  the  one  in  Fig.  9,  with  the  difference  of  its  ridge  width  being  increased  to  15 
pm.  The  spectra  do  not  broaden  much  even  at  very  high  pumping  current,  and  do  not 
show  the  Rabi  splitting.  As  discussed  in  Appendix  B,  the  effect  of  Kerr-lensing 
decreases  significantly  when  the  active  region  width  is  increased.  Thus,  this  is  a  strong 
evidence  that  the  saturable  absorption  needed  for  lowering  the  RNGH  threshold  is 
provided  by  the  Kerr-lensing  mechanism. 


7.  Experimental  results:  Buried  heterostructure  lasers 


The  second  type  of  lasers  we  have  tested  are  the  so-called  buried  heterostructure  lasers, 
in  which  an  insulating  Fe-doped  InP  layer  is  regrown  after  etching  of  the  ridges.  A  thick 
layer  of  electrically  plated  gold  is  deposited  on  top  of  the  ridges  as  top  contact  layer 
after  the  InP  regrowth  .  Figure  14  shows  the  cross  section  of  such  a  laser. 


14.  SEM  image  of  the  cross  section  of  a  buried  heterostructure  QC  laser. 


The  active  region  of  the  buried  heterostructure  lasers  tested  is  based  on  the 
four-quantum-well  design,  which  relies  on  a  double  phonon  resonance  to  achieve 
population  inversion  [39],  Fig.  15a  shows  the  V-I  and  L-I  characteristics  of  a  laser  with 
active  region  width  of  3  pm  and  wavelength  8.38  pm  operated  in  CW  at  room 
temperature,  and  b  shows  its  optical  spectra.  The  spectra  show  Rabi  splitting  similar  to 
the  ridge  laser  case,  indicating  the  RNGH  instability  in  this  narrow  buried 
heterostructure  laser.  The  spectral  splitting  and  twice  the  Rabi  frequency  are  plotted 
against  the  square  root  of  the  collected  output  power  in  Fig.  15  b.  A  good  agreement  is 
found  between  the  experimental  splitting  and  twice  the  calculated  Rabi  frequency. 


15.  (a)  V-I  and  LI  curves  and  (b)  optical  spectra  vs.  pumping  ratio  (j/jth)  above 
threshold  obtained  in  CW  at  300K  of  a  3  pm  wide  ridge  laser  emitting  at  8pm.  (c) 
Spectral  splitting  and  twice  the  Rabi  frequency  vs.  square  root  of  output  power 
collected  from  a  single  facet.  The  different  quantities  reported  on  the  graph  were 
deduced  from  the  experimental  data  show  in  (b).  The  dashed  line  is  the 
least-square  linear  fit  of  the  data. 


The  second-order  autocorrelation  trace  of  the  device  when  it  is  pumped  five  times 
above  threshold  at  80K  is  shown  in  Fig.  16.  The  ratio  between  the  maximum  of  the 
interference  fringes  and  the  background  is  close  to  8  to  3,  similar  to  the  ridge  laser  case. 
In  addition,  there  are  smaller  interference  fringes  within  one  cavity  roundtrip  period. 
Both  features  indicate  no  stable  pulsation  from  the  laser  (see  Appendix  A).  The 
microwave  spectrum  of  the  laser  output  shows  a  steady  peak  with  FWHM  of  700  kHz 
at  15.018  GHz.  This  is  shown  in  the  inset  of  Fig.  16. 
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16.  A  second-order  autocorrelation  trace  of  a  8  pm  wavelength  buried  heterostructure  QC  laser 
(wafer  #  3251)  under  the  condition  of  RNGH  instability.  (Inset)  Microwave  spectrum  of  the 
photocurrent  generated  by  the  same  laser  (measured  with  a  68  kHz  resolution  bandwidth). 


As  in  the  ridge  laser  case,  the  lowering  of  the  threshold  of  RNGH  instability  can  also  be 
attributed  to  the  Kerr-lens  type  saturable  absorption.  In  order  to  better  support  this  idea 
experimentally,  the  spectra  from  another  device  processed  from  the  same  wafer  as  in 
Fig.  15  but  with  a  wider  active  region  (7.5  pm)  are  measured.  Two-dimensional 
waveguide  simulations  indicate  a  much  weaker  Kerr-lensing  effect  in  these  QCLs,  due 
to  the  much  larger  ratio  of  active  region  width  to  wavelength  (see  Appendix  B).  The 
measured  optical  spectra  obtained  at  300K  in  CW  mode  are  shown  in  Fig.  17  a.  The 
envelopes  of  the  spectra  consist  of  multiple  peaks  whose  separation  is  independent  of 
the  pumping  current.  The  spectral  signatures  qualitatively  agree  with  the  numerical 
simulations  in  a  Fabry-Perot  cavity  without  a  saturable  absorber  (y  =0)  (Fig.  6). 
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17.  Optical  spectra  vs.  pumping  ratio  (j/jth)  above  threshold  obtained  in  CW  at  300K 
with  (a)  a  7.5  pm  wide  and  (b)  a  7.5  pm  buried  heterostructure  lasers  emitting  at  8pm. 


Further  increasing  the  active  region  width  tends  to  suppress  both  the  instabilities  caused 
by  RNGH  and  spatial  hole  burning.  Fig.  17  b  shows  the  spectra  of  a  device  also 
processed  from  the  same  wafer  but  with  an  even  wider  active  region  (10  pm)  at  77K. 
The  spectra  do  not  broaden  much  even  at  very  high  pumping  currents,  as  in  the  ridge 
laser  case.  Since  both  RNGH  instability  and  SHB  stem  from  nonlinear  effects,  they  are 
suppressed  when  the  intensity  of  the  field  in  the  cavity  is  lower.  Moreover,  when  the 
active  region  width  increases,  higher  transverse  modes  are  excited.  Different  transverse 
modes  have  different  propagation  constants  p,  and  thus  form  different  gain  gratings 
which  tend  to  wash  out  the  effect  of  spatial  hole  burning. 


8.  Temperature  effects 


18.  Optical  spectra  vs.  pumping  ratio  (j/jth)  above  threshold  obtained  in  CW  of  a 
buried  heterostructure  emitting  at  8.38pm  at  different  temperatures:  (a)  80K,  (b) 

150K,  (c)  200K  and  (d)  300K. 


Temperature  also  plays  an  important  role  in  the  shape  of  the  spectrum.  To  illustrate  this 
point,  we  now  present  the  spectra  at  different  temperatures  for  the  >,=8.38  pm  buried 
heterostructure  laser  with  3  pm  active  region  width  (the  same  as  in  Fig.  15). 

From  the  spectra  in  Fig.  18,  it  is  clear  that  at  lower  temperatures  the  spectra  are 
dominated  by  spatial  hole  burning,  showing  multiple  peaks  independent  of  pumping 
and  no  significant  Rabi  splitting.  As  the  temperature  increases,  the  Rabi  splitting 
becomes  more  evident  and  finally  the  RNGH  instability  takes  over. 

Intuitively,  the  effect  of  the  temperature  on  the  nature  of  the  multimode  regime  in  QCLs 
seems  to  come  from  carrier  diffusion.  The  lifetime  of  the  gain  grating  Tg  is  defined  as 
Tg'  =  T~'  +4  k2D .  The  diffusion  coefficient  D  of  the  gain  grating  is  proportional  to  the 

temperature,/)  =  jukBT / q  ,  where  p  here  is  the  carrier  mobility,  kB  is  the  Boltzmann’s 
constant,  T  is  the  temperature  and  q  is  the  carrier  charge.  Therefore  at  higher 
temperature  carrier  diffusion  would  reduce  spatial  hole  burning  and  thus  reveal  the 
RNGH  instability.  However,  with  a  mobility  of  16000  cm2/sec/V  and  k=2.25*104  cm'1 
(which  corresonds  to  a  vacuum  wavelength  of  8.38pm),  4k2D  ~  0.2  THz  at  77K  and 
4k2D  ~  0.8  THz  at  300K,  both  significantly  smaller  than  Tf1  ~  0.6  THz  at  77K  and  Tf1 
~  2  THz  at  300K.  Thus,  carrier  diffusion  is  unlikely  to  be  the  reason  for  the  temperature 
dependence.  Although  the  temperature  effect  is  not  entirely  understood  to  us  at  this 
point,  one  possiblity  is  temperature-dependent  saturable  absorption.  The  QCL  injector 
consists  of  many  energy  levels  which  can  be  thermally  populated.  It  is  not  surprising 
that  if  any  two  higher  levels  in  the  injector  are  closely  resonant  with  the  laser  transition, 
and  that  they  will  form  a  two-level  saturable  absorber. Therefore  in  this  case  saturable 
absorption  is  stronger  at  higher  temperature  and  makes  the  RNGH  more  easily 
observable. 


9.  Conclusion 


This  paper  provides  a  thorough  account  of  different  multimode  regimes  in  QCLs.  It  was 
found  that  two  key  mechanisms  which  govern  the  multimode  regimes  in  QCLs  are  a 
coherent  instability  similar  to  the  RNGH  instability,  and  spatial  hole  burning.  Both 
mechanisms  are  enhanced  due  to  the  large  dipole  moment  p  of  the  laser  transition, 
which  results  in  the  unusually  fast  gain  recovery  in  QCLs:  The  RNGH  instability  is 
enhanced  because  the  Rabi  frequency  scales  as  p,  and  therefore  the  Rabi  splitting  can 
be  resolved  by  the  comb  of  modes  supported  by  the  cavity.  SHB  is  enhanced  because 
carrier  diffusion  is  slower  than  the  gain  recovery,  and  thus  leaves  the  gain  grating 
intact. 

Due  to  the  fast  gain  recovery,  conventional  mode  locking,  with  one  pulse  per  roundtrip, 
is  suppressed.  In  order  to  achieve  conventional  mode  locking  in  QCLs,  one  needs  to 
design  a  QCL  with  a  slower  gain  recovery,  such  that  Ti  becomes  longer  than  or 
comparable  to  the  cavity  roundtrip.  Efforts  in  this  direction  are  currently  underway. 
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Appendix  A.  Interferometric  autocorrelation 

If  we  write  the  electric  field  as  function  of  time  as  E(t)elC0t  ,  the  two-photon 
interferometric  autocorrelation  is  given  by 
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We  define  the  background  as 


From  Eq.  (65)  one  can  see  that  I(0)=8Ib.  We  now  assume  that  E(t)  is  an  isolated  pulse, 
and  for  the  sake  of  simplicity  we  assume  that  E(t)  is  nonzero  only  over  an  interval  T. 
Then  for  x  >T  all  terms  in  Eq.  (65)  which  include  both  t  and  t+x  vanish.  Therefore 
I(x>T)=Ib.  In  other  words,  for  a  pulse  one  has 
1(0) 
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Now  let  us  assume  that  E(t)  is  a  complex  stationary  random  process  whose  phase  at 
each  point  in  time  is  uniformly  distributed  over  the  interval  [0,  2n].  Let  us  assume  for 
simplicity  that  there  is  T  such  that  for  x>T,  E(t)  and  E(t+x)  are  statistically  independent. 
One  can  see  that  the  mean  value  of  the  third,  the  fourth  and  the  fifth  line  of  Eq.  (65) 
vanish.  The  second  line  of  Eq.  (65)  equals  31b-  We  therefore  obtain  that  for  a  random 
process 
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The  ratio  between  the  peak  and  the  background  in  the  autocorrelation  function  is 
therefore  8:1  for  an  isolated  pulse,  that  is  when  the  phases  of  all  modes  are  all  zero,  and 
8:3  for  modes  with  completely  random  phases. 


19.  Simulated  interferometric  autocorrelation,  corresponding  to  the  same  parameters  as  in  Fig.  7. 

The  value  of  p  is  denoted  at  each  plot. 

Fig.  19  shows  simulated  interferometric  autocorrelation  traces  for  the  same  parameters 
as  in  Fig.  6.  For  p=l  the  phases  of  the  modes  are  random,  and  the  peak  to  background 
ratio  is  8:3.  When  the  saturable  absorber  is  absent,  y=0,  the  ratio  in  the  autocorrelation 
is  also  8:3.  When  the  pumping  is  higher  (Fig.  19  below),  the  saturable  absorber  induces 
some  phase  relationships  between  the  modes,  and  the  peak  to  background  ratio  is  about 
8:2.  Note  that  in  the  autocorrelation  in  Fig.  12,  the  ratio  is  slightly  greater  than  8:3.  The 
autocorrelation  also  has  a  nontrivial  structure,  with  a  peak  at  half  the  cavity  roundtrip 
time.  We  were  not  able  to  reproduce  this  structure  in  the  simulations. 


Appendix  B.  Kerr  lensing  in  QCLs 

This  section  elaborates  on  the  analysis  of  the  Kerr  lensing  effect  (including  soft  Kerr 
lensing)  in  QCLs.  The  mechanism  of  Kerr  lensing  can  be  understood  as  follows:  above 
lasing  threshold,  there  is  net  gain  in  the  active  region  which  compensates  the  mirror 
losses.  Assuming  a  nonlinear  refractive  index  A nNL  =  n2  ■  I  in  the  active  region,  it 

causes  self  focusing  of  the  transverse  mode,  which  results  in  a  net  increase  in  the  modal 
gain  due  to  an  increased  overlap  with  the  active  region  and  a  decreased  overlap  with 
lossy  waveguide  cladding.  A  stronger  intensity  leads  to  stronger  Kerr  lensing,  thus 
forming  an  intensity-dependent  saturable  absorber. 

The  active  region  width  is  expected  to  play  an  important  role  in  the  Kerr  lensing  effect. 
The  narrower  the  active  region,  the  less  confined  is  the  transverse  mode,  which  results 
in  a  bigger  increase  of  the  modal  gain  due  to  Kerr  lensing.  Therefore  the  saturable 
absorber  coefficient  y  is  expected  to  be  larger  in  lasers  with  narrower  active  region. 
We  performed  FDTD  waveguide  simulations  using  the  commercial  software 
BeamPROP.  We  simulated  buried  heterostructure  QCLs  (wafer  #  3251)  with  various 
active  region  widths.  The  software  allows  us  to  assign  complex  refractive  indices  to 
each  layer,  and  the  effective  modal  index  is  calculated  by  the  software.  The  detailed 
steps  of  our  Kerr  lensing  simulation  are  as  follows: 

1 .  We  first  assign  a  net  gain  (i.e.,  a  negative  imaginary  part  of  nactive)  in  the 
active  region,  making  the  total  modal  gain  equal  to  the  total  mirror  loss.  In  our 
case,  given  the  reflectivity  r=0.27  at  our  laser  facets,  a  3mm  long  laser  hasa  total 
mirror  loss  of  4  cm'1. 

2.  We  assign  different  values  of  A nNL  =  n2  ■  I  to  the  real  part  of  nactive- 
After  running  the  simulation,  we  get  the  net  modal  gain  for  each  nNL- 

3.  We  plot  the  modal  gain  versus  A nNL  .  For  example,  Fig.  20  shows  the 
plots  for  3  pm  wide  and  5  pm  wide  active  regions. 

4.  In  order  to  relate  to  our  theoretical  model,  the  change  in  modal  gain  is 
attributed  entirely  to  change  in  modal  losses. 

5.  The  change  in  modal  loss  equals  to  y |2sj 2 ,  which  is  also  proportional  to 

the  intensity  I  Thus,  from  the  slope  of  the  plot  of  modal  gain  vs. 
A nNL  =  n2- 1  ,we  can  extract  the  ratio  of  y/n2. 

6.  For  a  given  active  region  width,  we  choose  the  value  of  y  that  brings  the 
threshold  of  RNGH  instability  down  close  to  our  experimental  data.  Then  we 
can  determine  the  value  of  n2  from  the  ratio  of  y/n2. 

7.  With  the  fixed  value  of  n2,  we  repeat  the  waveguide  simulation  with 
different  active  region  widths  and  plot  y  vs.  the  active  region  width. 


20.  Simulation  results  of  the  modal  gain  vs.  the  nonlinear  refractive  index  [Trial  mode]  for  a 
3  Jim  wide  and  a  5  Jim  wide  active  region  in  wafer  #  3251. 

For  an  active  region  width  of  3  pm,  the  RNGH  threshold  is  dramatically  lowered  from 
nine  times  the  lasing  threshold  to  less  than  twice  the  lasing  threshold  (see  Fig.  11). 
From  Fig.  2,  the  required  y  is  about  2*  10"9  cm/V2,  from  which  we  obtain  that  n2  is  about 
2*10'  cm  /W.  The  estimated  y  vs.  active  region  width  is  plotted  in  Fig.  21.  As  the 
active  region  width  increases,  y  decreases,  and  when  the  active  region  is  as  wide  as  7.5 
pm,  y  decreases  to  1*1  O’9  cm/V2.  This  pushes  the  RNGH  instability  threshold  to  above 
twice  the  lasing  threhsold.  This  may  explain  why  we  see  only  instabilities  resulting 
from  spatial  hole  burning  in  wider  buried  heterostructure  lasers  (see  Fig.  17). 


21.  Simulation  results  of  the  saturable  absorber  coefficient  y  vs.  the  active  region  width  in  wafer  # 

3251  assuming  n2=2*10'8  cm2/W. 

For  ridge  lasers,  the  lossy  gold  contact  on  the  sidewalls  of  the  laser  ridges  will  make  the 
saturable  absorber  even  stronger.  This  may  explain  the  absence  of  spatial  hole  burning 
in  ridge  lasers. 

There  are  several  possibilities  which  might  explain  the  origin  of  n2  in  the  active  region. 
First,  the  material  we  use  has  a  thermal  index  change  of 

An  /AT  =  1.55x10_4AT"1,  and  when  the  laser  is  operating  in  CW,  the  active  region 
usually  heats  up  by  tens  of  degrees.  Second,  the  complex  states  in  the  injector  can  all 
contribute  to  both  real  and  imaginary  parts  of  x<3)  at  the  laser  wavelength,  which  results 
in  both  n2  and  y.  The  latter  effect  should  be  enhanced  with  increasing  temperature  since 
electrons  populate  a  larger  amount  of  absorbing  states  at  higher  temperatures.  Third, 
although  not  very  significant,  there  is  always  bulk  n2  of  the  material. 
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We  present  the  pulsed  operation  at  room  temperature  of  different  strained  InGaAs/AlInAs 
quantum-cascade  lasers  grown  by  low-pressure  metalorganic  vapor-phase  epitaxy.  Devices  based  on 
a  bound-to-continuum  transition  design  have  threshold  current  densities  in  pulsed  mode  as  low  as 
1.84  kA/cm2  at  300  K.  Identical  lasers  grown  at  higher  rate  (0.5  nm/s)  also  have  threshold  current 
densities  lower  than  2  kA/cm2  at  300  K.  Buried  hetero structure  lasers  based  on  a  double  phonon 
resonance  design  were  operated  in  continuous  mode  up  to  280  K.  Overall,  the  performance  obtained 
from  strained  quantum  cascade  lasers  deposited  by  metalorganic  vapor-phase  epitaxy  are 
comparable  with  that  of  similar  structures  grown  by  molecular  beam  epitaxy.  ©  2006  American 
Institute  of  Physics.  [DOI:  10.1063/1.2166206] 


Quantum  cascade  lasers  (QCLs)  are  the  light  source  of 
choice  for  many  potential  applications  such  as  gas  sensing  in 
the  midinfrared.1  State-of-the-art  QCLs  have  indeed  demon¬ 
strated  large  tunability,  low-threshold,  and  high-power  in 
continuous  mode  (cw)  at  room  temperature.2,3  The  devices 
which  exhibit  the  best  performance  have  been  so  far  grown 
by  molecular  beam  epitaxy  (MBE).  This  fabrication  tech¬ 
nique  provides  several  advantages,  such  as  the  ease  to  form 
abrupt  interfaces,  but  is  not  suitable  for  low-cost,  high 
throughput  production.  On  the  contrary,  low-pressure  metal¬ 
organic  vapor-phase  epitaxy  (MOVPE)  is  a  well-established 
technology  fitting  the  requirements  of  large  scale  production, 
especially  since  high  growth  rates  are  possible.4-6  However, 
it  is  still  necessary  to  demonstrate  that  MOVPE-  and  MBE- 
grown  QCLs  can  have  similar  performance  particularly  since 
the  control  of  interface  abruptness  is  more  challenging  by 
MOVPE.  QCLs  fabricated  using  the  latter  technique  have 
already  been  operated  recently  at  room-temperature  in  pulsed 
mode.  The  threshold  current  densities  for  these  lattice- 
matched  devices  were  low,  i.e.,  between  2.4  and  3  kA/cm2  at 
300  K  and  the  emission  wavelength  longer  than  7.2  /zm.4,6 
Structures  at  shorter  wavelengths  require  in  order  to  achieve 
high  performance  at  room  temperature,  the  use  of  strain- 
balanced  material  because  of  the  large  band  discontinuity 
available.7 

This  letter  reports  on  the  fabrication  of  strained  QC  la¬ 
sers  grown  by  MOVPE  working  in  pulsed  mode  above 
320  K.  Threshold  current  density  as  low  as  1.84  kA/cm2  at 
300  K  and  characteristic  temperature  T0  close  to  200  K  were 
obtained  from  devices  based  on  a  bound-to-continuum  de¬ 
sign.  Another  sample  based  on  a  double  phonon  resonance 
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design  worked  in  cw  up  to  280  K.  More  than  80  mW  output 
power  was  delivered  by  the  device  in  cw  mode  at  243  K. 

Three  strained  QC  structures  were  deposited  on  highly 
doped  substrate  by  low  pressure  (76  Torr)  MOVPE.  The  pa¬ 
rameters  were  identical  to  those  used  for  the  growth  of  the 
lattice-matched  QCLs  reported  in  Refs.  5  and  6.  The  flows  of 
triethylgallium  and  trimethylaluminium  were  simply  ad¬ 
justed  to  form  the  strained  In06Ga04As  and  Al056In044As 
alloys.  The  active  region  of  the  first  two  samples  investigated 
is  based  on  a  bound-to-continuum  design  previously  reported 
by  Blaser  et  al?  and  was  deposited  at  a  slow  (0.1  nm/s)  and 
a  fast  (0.5  nm/s)  rate.  The  doping  in  the  injector  was  kept 
low  (lXl017cm3)  to  limit  the  waveguide  losses.  Thirty 
stages  were  grown  and  embedded  between  0.3-/zm-thick  In- 
GaAs  guiding  layers  doped  3  X  1016  cm3,  in  addition  to 
3-/zm-thick  InP  cladding  layers  doped  1  X  1017  cm3.  Step- 
graded  30-nm-thick  InGaAsP  layers  were  inserted  at  the 
InP-InGaAs  interfaces  to  reduce  series  resistance.  The 
growth  ended  with  a  0.5 -yam- thick  InP  layer  doped 
1  X  1019  cm3  for  plasmon-enhanced  optical  confinement  and 
reduced  losses,  followed  by  highly  doped  contact  layers.  The 
deposition  rate,  except  in  the  active  region  was  kept  constant 
at  a  value  of  0.5  nm/s.  The  third  structure  studied  in  this 
letter  replicates  the  lasers  described  in  Refs.  8  and  9  based  on 
a  double  phonon  resonance  design.  Our  sample  comprises  30 
stages  and  the  waveguide  is  similar  to  the  one  described 
earlier. 

Figure  1  shows  the  x-ray  diffraction  spectra  of  the  two 
QCLs  based  on  a  bound-to-continuum  design,  together  with 
the  results  of  simulations.  The  comparison  between  the  dif¬ 
ferent  patterns  clearly  shows  very  similar  diffraction  spectra, 
no  noticeable  degradation  of  the  crystal  quality  as  the  growth 
rate  changes  and  an  excellent  reproducibility  of  the  structure 
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FIG.  1 .  Comparison  between  simulations  and  x-ray  diffraction  spectra  of  the 
fast-  and  slow-grown  bound-to-continuum  QCL  measured  along  the  (004) 
crystal  direction.  The  curves  are  shifted  for  clarity. 

from  one  run  to  the  other.  For  both  samples,  the  zeroth  order 
peak  due  to  the  periodicity  of  the  superlattice  is  very  close  to 
the  InP  substrate  peak,  which  indicates  only  a  slight  lattice 
mismatch  in  the  structures.  The  average  length  of  the  QC 
stages  can  be  deduced  from  the  different  spectra.  The  period 
of  the  active  region  is  49.3  and  49.7  nm  for  the  slow-  and 
fast-grown  structures,  respectively,  which  is  in  both  cases 
thinner  than  the  targeted  thickness  (50.6  nm).  This  observa¬ 
tion  is  consistent  with  the  fact  that  the  emission  wavelength 
of  our  lasers  is  blueshifted  with  respect  to  the  original  design 
[4.98  jum  vs  5.4  /z m,  see  the  inset  of  Fig.  2(a)]. 

Ridge  waveguide  lasers  were  fabricated  using  conven¬ 
tional  processing  techniques.  The  samples  were  cleaved 
manually,  soldered  ridge  side-up  with  indium  onto  gold- 
plated  copper  heatsinks  and  finally  wire  bonded.  After  fabri¬ 
cation,  the  lasers  were  mounted  into  a  temperature-controlled 
LN2-cooled  flow  cryostat.  Current  pulses  (163  ns)  at  a  fre¬ 
quency  of  92  kHz,  corresponding  to  a  duty  cycle  of  1.5%, 
were  supplied  to  the  lasers.  For  measurements  in  continuous 
mode,  the  devices  were  mounted  on  a  Peltier  cooler  installed 
in  a  purged  box.  The  light  from  the  laser  facet  was  collected 
by  // 1  optics  and  sent  either  onto  a  calibrated  thermopile 
detector  for  power  measurements  or  into  a  Fourier-transform 
infrared  spectrometer  equipped  with  a  deuterated  triglycine 
sulphate  detector  for  spectral  characterization. 

Figure  2(a)  shows  the  voltage  and  the  light-intensity  ver¬ 
sus  current  (VI  and  LI)  curves  of  a  1.9-mm-long  and 
16-/zm-wide  laser  fabricated  from  the  slow-grown  material. 
The  maximum  emitted  peak  power  was  156  and  51  mW  with 
a  slope  efficiency  dPIdl  equal  to  564  and  380  mW/A  at, 
respectively,  243  and  300  K.  The  inset  of  Fig.  2(a)  shows 
typical  spectra  obtained  close  to  threshold. 

The  threshold  current  density  yth  measured  in  pulsed 
mode  is  plotted  as  function  of  temperature  in  Fig.  2(b)  for 
the  fast-  and  slow-grown  samples.  In  both  cases,  the  lasers 
worked  at  temperatures  above  320  K  and  the  characteristic 
temperature  T0,  which  describes  the  behavior  of  the  thresh¬ 
old  current  versus  temperature  was  close  to  200  K.  For  the 
slow-grown  device,  jth  ranges  from  1.39  to  2.07  kA/cm2  be¬ 
tween  243  and  320  K  and  is  close  to  the  values  obtained  in 
Ref.  3.  In  this  letter  indeed,  a  threshold  current  density  of 
1.38  kA/cm2  was  reported  for  a  3-mm-long  device  operated 
in  continuous  mode  at  a  heatsink  temperature  Tsink  of  243  K. 
The  latter  corresponds  to  a  temperature  in  the  active  region 


(a)  Current  [A] 


FIG.  2.  Summary  of  the  data  obtained  with  the  lasers  based  on  a  bound-to- 
continuum  transition  design  in  pulsed  mode,  (a)  VI /LI  curves  and  optical 
spectra  close  to  threshold  (inset)  taken  at  243  and  300  K  with  the  slow 
grown  sample,  (b)  Threshold  current  density  vs  temperature  for  the  slow- 
and  fast-grown  QCLs.  Note  that  the  latter  has  a  cavity  length  31%  longer 
than  the  slow-grown  laser.  The  curves  correspond  to  the  usual  exponential  fit 
jth=jo  exp (77 r0). 

Pelec  of  3.75  W  and  the  thermal  resistance  of  Rth  6.3  K/W 
measured  by  B laser  et  al.  Assuming  that  heating  effects  are 
negligible  in  pulsed  mode,  we  can  expect  for  their  device  a 
pulsed  threshold  of  1.38  kA/cm2  at  about  266  K.  This  is 
only  15%  lower  than  the  threshold  current  density  of  the 
slow-grown  QCL  investigated  in  the  present  letter.  This  dif¬ 
ference  would  be  even  smaller  lasers  of  equal  length  were 
compared. 

As  shown  in  Fig.  2(b),  the  threshold  current  density  of 
the  fast-grown  lasers  is  less  than  2  kA/cm2  up  to  304  K  and 
the  characteristic  temperature  T{)  reaches  185  K.  This  result 
is  remarkable  given  the  high  deposition  rate  (0.5  nm/s)  used 
during  the  growth  of  the  sample.  However,  the  difference  in 
terms  of  performance  between  the  slow-  and  fast-grown 
QCLs  is  not  negligible,  although  the  curves  in  Fig.  2(b)  sug¬ 
gest  the  contrary.  It  is  due  to  the  fact  that  the  cavity  length  of 
the  fast-grown  QCL  was  longer  (2.52  mm)  than  in  the  case 
of  the  slow-grown  sample.  If  both  lasers  would  have  been 
cleaved  in  exactly  2-mm-long  bars,  the  difference  in  thresh¬ 
old  current  density  would  be  about  20%  higher.  Note  that  the 
measured  slope  efficiency  is  only  232  mW/A  for  the  fast 
grown  samples. 

The  QC  sample  based  on  a  double  phonon  resonance 
design  was  processed  into  narrow-stripe  buried  heterostruc¬ 
ture  lasers.  In  this  case,  the  active  region  is  surrounded  by 
insulating  InP  doped  with  Fe,  lowering  the  thermal  resistance 
of  the  device  and  allowing  continuous  wave  operation  at 


^act-^sink+^eiec‘^th  dose  to  266  K,  given  the  electric  power  high  temperatures.  To  reduce  the  thermal  resistance  even  fur- 
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(b)  Temperature  [K] 

FIG.  3.  Summary  of  the  data  obtained  with  the  structure  based  on  a  four 
quantum  wells  transition  design  in  pulsed  and  continuous  mode,  (a)  VI  and 
LI  taken  at  different  temperatures  in  continuous  mode  after  the  evaporation 
of  a  HR  coating  on  the  backfacet.  The  inset  shows  the  optical  spectra  mea¬ 
sured  close  to  threshold  at  243  K.  (b)  Threshold  current  density  as  a  func¬ 
tion  of  the  heatsink  temperature  in  pulsed  and  continuous  mode,  respectively 
without  and  with  a  high-reflection  coating.  The  curves  correspond  to  the 
usual  exponential  fit  y'th=y0  exp(77r0). 

ther,  a  thick  layer  of  gold  was  electroplated  onto  the  usual 
Ti/Au  top  contact.  The  best  results  were  obtained  with  a 
7-/zm-wide  and  2-mm-long  laser  mounted  junction-up, 
which  operated  at  temperatures  above  320  K  in  pulsed  mode. 
The  measured  slope  efficiency  with  a  1.5%  duty  cycle 
reached  636  mW/A  at  243  K  and  339  mW/A  at  320  K.  The 
device  could  be  operated  in  continuous  mode  at  243  K.  The 
temperature  range  accessible  was,  however,  limited  since  the 
cw  threshold  current  density  (3.29  kA/cm2  at  243  K)  was 
fairly  close  to  the  injection  current  at  which  the  roll-over  of 
the  output  power  takes  place  (3.74  kA/cm2  at  243  K).  The 
maximum  intensity  was  6.6  mW  and  the  slope  efficiency 
172  mW/A.  A  high-reflection  (HR)  coating  consisting 
of  175  nm  A1203  and  30  nm  Au  was  deposited  on  the 
backfacet  of  the  laser  and  helped  to  reduce  the  cw  threshold 
current  density.  The  value  of  the  latter  dropped  to 
2.61  kA/cm2  at  243  K  and  the  maximum  cw  temperature 
operation  of  the  laser  increased  up  to  280  K.  Figure 
3(a)  shows  the  cw  VI-LI  curves  of  the  device  after  the 
evaporation  of  the  HR  coating.  The  intensity  reached  81  mW 
(dP/dI=651  mW/A)  at  243  K  and  was  still  8  mW  ( dP/dl 
=  251  mW/A)  at  273  K.  The  inset  of  Fig.  3(a)  shows  a  typi¬ 
cal  spectrum  obtained  if  the  laser  is  operated  slightly  above 
threshold  in  cw  mode.  At  243  K,  the  emission  wavelength  is 
1954  cm-1  (5.11  /mi),  which  is  in  quite  good  agreement 
with  the  expected  value  of  (5.3  /mi).8 

The  threshold  current  densities  obtained  in  pulsed  and 
cw  mode  at  different  temperatures  are  displayed  in  Fig.  3(b). 
Note  that  the  measurements  shown  were  taken  without 
(pulsed  mode)  and  with  (cw)  the  HR  coating  on  the  back- 


facet  of  the  laser.  From  an  exponential  fit  of  the  data,  the 
characteristic  temperature  T{)  can  be  deduced.  It  yields  207 
and  111  K  in  pulsed  and  cw  mode,  respectively.  The  thermal 
resistance  Rth  of  the  device  can  be  estimated  at  243  K  by 
comparing  the  threshold  current  density  in  cw  obtained  with¬ 
out  HR  coating  (3.29  kA/cm2)  and  the  pulsed  measurements 
displayed  in  Fig.  3(b).  This  leads  to  a  value  of  Rth  equal  to 
14.5  K/W  corresponding  to  a  thermal  conductance  Gth 
=  516  W/K  cm-2.  From  the  change  of  current  density  before 
and  after  the  deposition  of  the  HR  coating,  the  waveguide 
losses  can  be  estimated  and  yield  9.7  cm-1.  The  perfor¬ 
mances  reported  in  Figs.  3(a)  and  3(b)  are  similar  to  those 
reported  in  Refs.  8  and  9  although  a  direct  comparison  is  not 
as  straight-forward  as  in  the  case  of  the  bound-to-continuum 
samples.  Slightly  lower  performance  is  however,  anticipated 
for  the  MOVPE-grown  lasers  since  the  linewidth  of  the  elec¬ 
troluminescence  peak  measured  at  room  temperature  is 
broader  (36  meV)  than  the  corresponding  value  (25  meV) 
found  in  Ref.  8  and  9. 

In  conclusion  our  results  demonstrate  the  realization  of 
strained  QCLs  based  on  a  bound-to-continuum  design  which 
operate  in  pulsed  mode  above  320  K  with  a  very  low  thresh¬ 
old  current  density.  Devices  based  on  a  double  phonon  reso¬ 
nance  design  and  processed  into  buried  heterostructures  la¬ 
sers  worked  in  continuous  wave  up  to  280  K.  The  results 
presented  in  this  letter  are  comparable  with  those  of  equiva¬ 
lent  structures  grown  by  molecular  beam  epitaxy.  The  per¬ 
formance  level  can  be  improved  in  many  ways,  including  the 
optimization  of  the  processing  of  buried  heterostructure  de¬ 
vices  as  well  as  the  doping  in  the  active  region.  It  will  allow 
without  doubts  the  cw  operation  at  300  K  of  these  lasers. 
Another  remarkable  finding  is  the  low  threshold  current  den¬ 
sities  (less  than  2  kA/cm2)  measured  at  room  temperature 
with  lasers  deposited  at  a  high  growth  rate  of  0.5  nm/s. 
Altogether,  the  results  presented  in  this  letter  demonstrate 
that  MOVPE  is  a  viable  solution  to  grow  high  performance 
low-threshold  strained  QCLs  operating  in  cw  at  room 
temperature. 
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High-power  quantum  cascade  lasers  (QCLs)  working  in  continuous  wave  (cw)  above  400  K  are 
presented.  The  material  was  grown  by  low-pressure  metal  organic  vapor-phase  epitaxy  and 
processed  into  narrow  buried  hetero structure  lasers.  A  cw  output  power  of  204  mW  was  obtained  at 
300  K  with  an  8.38  jam  wavelength,  3  mm  long  and  7.5  jam  wide  coated  laser.  The  device  operates 
in  cw  mode  above  400  K,  which  exceeds  the  previous  maximum  cw  temperature  operation  of  QCLs 
by  approximately  60  K.  Preliminary  reliability  data  obtained  by  accelerated  aging  tests  indicate  a 
remarkable  robustness  of  the  lasers.  ©  2006  American  Institute  of  Physics. 
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It  has  been  shown  in  the  last  few  years  that  quantum 
cascade  lasers  (QCLs)  are  an  ideal  light  source  for  many 
chemical  sensing  applications  in  the  midinfrared.1,2  Indeed 
stable  single  mode  operation,  which  is  highly  desirable  for 
these  applications  can  be  achieved  with  QCLs  by  means  of 
an  external  cavity3  or  by  fabricating  distributed  feedback 
lasers,4  for  example.  In  addition,  QCLs  can  work  in  continu¬ 
ous  wave  (cw)  at  temperatures  above  300  K  with  high  output 
power,  as  discussed  in  Refs.  4  and  5.  A  recent  article  reported 
a  cw  power  of  446  mW  at  292  K  and  a  maximum  tempera¬ 
ture  operation  of  333  K  for  a  15  jam  wide,  3  mm  long  device 
emitting  at  6  /zm.6  As  the  commercial  significance  of  QCLs 
increases,  the  production  of  reliable  and  potentially  low-cost 
devices  becomes  increasingly  important.  The  question  of 
whether  or  not  high  performance  can  be  obtained  from  QCLs 
grown  by  metal  organic  vapor-phase  epitaxy  (MOVPE)  is 
central  since  this  technique  is  a  widely  established  platform 
for  the  high- volume  production  of  reliable  semiconductor  la¬ 
sers.  This  question  was  partially  answered  recently  by  the 
demonstration  of  room  temperature  cw  operation  with  lattice 
matched  QCLs  grown  by  MOVPE.7  The  devices  tested, 
which  were  based  on  a  basic  three  quantum  well  active  re¬ 
gion  design,  worked  cw  up  to  320  K  with  an  output  power  of 
more  than  20  mW  at  300  K. 

In  this  letter  we  report  the  fabrication  of  MOVPE-grown 
QCLs  having  performance  comparable  to  the  best  results  ob¬ 
tained  with  MBE  material.6  The  maximum  collected  output 
power  from  a  7.5  jum  wide,  3  mm  long  QCL  is  as  high  as 
204  mW  at  300  K  after  the  evaporation  of  a  high-reflection 
(HR)  coating.  This  device  operated  in  continuous  mode  up  to 
400  K  with  more  than  10  mW  of  output  power.  Note  that 
throughout  this  letter,  the  values  of  the  output  power  and  the 
slope  efficiency  are  not  corrected  for  the  collection  efficiency 
of  our  setup,  which  is  estimated  to  be  70%. 

The  QCLs  investigated  here  were  grown  by  low- 
pressure  MOVPE  in  a  standard  reactor  equipped  with  a 
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purged  and  pressure-balanced  switching  manifold  and  a 
close-spaced  showerhead  injector.  More  details  about  the 
growth  can  be  found  in  Refs.  8  and  9.  The  active  region  of 
the  samples  studied  is  based  on  a  double-phonon  resonance 
design  identical,  except  for  the  doping  level  in  the  injector,  to 
the  design  reported  by  Hofstetter  et  al.10  35  states  were 
grown  at  a  slow  rate  of  0.1  nm/s.  The  doping  in  the  injector 
was  kept  low  (1  X  1017  cm3)  to  limit  the  optical  losses.  The 
waveguide  surrounding  the  active  region  consists  of  two  In- 
GaAs  layers  (thickness  of  0.52  jum,  Si:doping  of  3 
X  10 16  cm3)  and  two  thick  InP  cladding  layers  (thickness  of 
3.5  jam,  Si:doping  of  1  X  1017  cm3)  grown  at  a  fast  rate  of 
0.5  nm/s.  The  growth  sequence  ended  with  the  deposition  of 
a  plasmon-conhnement  layer  (InP,  thickness  of  0.5  jam, 
Si: doping  of  5  X  10 18  cm3)  and  20  nm  n+  InGaAs  capping 
layer. 

Buried  hetero  structure  lasers  were  processed  using 
conventional  fabrication  techniques.  Ridges  with  a  width  of 
either  3  or  7.5  jam  were  etched  by  a  combination  of  wet  and 
dry  etching.  Fe-doped  InP  was  then  regrown  to  planarize  the 
structure,  helping  to  decouple  laterally  the  optical  mode  from 
the  lossy  metal  contacts  and  to  lower  the  thermal  resistance 
of  the  device  while  minimizing  leakage  currents.  Electrical 
contacts  are  provided  by  Ti/Au  metallizations  evaporated 
directly  onto  the  n+  InGaAs  capping  layer  and  the  Fe:InP 
surrounding  the  laser  ridge.  An  additional  5  /zm  thick  gold 
layer  was  subsequently  electroplated  to  further  improve  heat 
dissipation.  The  devices  were  finally  mounted  ridge  side  up 
onto  copper  blocks  and  wire  bonded.  The  lasers  were  tested 
inside  a  Peltier-cooled  box  with  a  CaF2  window.  The  tem¬ 
perature  of  the  heat  sink  was  monitored  with  a  Si  diode 
placed  near  the  QCLs.  The  light  from  one  facet  was 
collected  by  // 1  optics  and  a  calibrated  thermopile  detector 
was  used  to  perform  power  measurements.  The  optical  spec¬ 
tra  were  recorded  with  a  Nicolet  860  Fourier  transform 
interferometer. 

Figure  1  shows  the  voltage  and  the  light  intensity  versus 
current  ( V-I  and  L-I )  curves  of  3  mm  long  lasers  having, 


0003-6951 /2006/88(20)/201 1 1 5/3/$23.00  88,  201115-1  ©  2006  American  Institute  of  Physics 

Downloaded  25  May  2006  to  128.103.60.225.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://apl.aip.org/apl/copyright.jsp 


201115-2 


Diehl  et  al. 


Appl.  Phys.  Lett.  88,  201115  (2006) 


Threshold  current  density  [kA/cm2] 

0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0 


(b) 


Current  [A] 

Threshold  current  density  [kA/cm2] 

0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0 


CD 

I 


FIG.  1.  cw  L-I  and  V-I  curves  of  HR  coated  3  mm  long  lasers  with  (a) 
3  jmm  and  (b)  7.5  jum  ridge  widths.  The  heat  sink  temperature  was  varied 
between  240  and  400  K  in  steps  of  20  K.  The  optical  power  data,  taken 
from  one  laser  facet,  are  not  corrected  for  the  estimated  70%  collection 
efficiency  of  our  setup. 


Temperature  [K] 

FIG.  2.  (a)  Threshold  current  density  as  a  function  of  heat  sink  temperature 
for  the  HR  coated,  3  mm  long  lasers  with  3  jum  (circles)  and  7.5  jum 
(squares)  ridge  width  in  pulsed  (dashed  line)  and  cw  modes  (solid  line).  The 
experimental  points  (open  symbols)  are  fitted  by  the  expression  jth=j0  eTIT°. 
(b)  Slope  efficiency  vs  temperature  deduced  from  the  L-I  curves  shown  in 
Fig.  1(b)  close  to  threshold.  Inset:  high-resolution  cw  spectra  as  a  function 
of  injection  current  measured  with  the  3  gm  wide  QCL  without  HR  coating 
at  a  constant  temperature  of  300  K. 


respectively,  3  and  7.5  jum  ridge  widths.  The  measurements 
were  performed  with  the  devices  operating  in  cw  mode  and 
at  temperatures  ranging  from  240  to  400  K.  A  HR  coating 
consisting  of  200  nm  A1203  and  30  nm  Au  was  evaporated 
on  the  backfacet  of  both  devices  and  nearly  doubled  the 
maximum  collected  output  power  Pmax  of  both  the  narrow 
and  the  broad  laser.  The  value  of  Pm ax,  measured  at  room 
temperature  with  the  7.5  jum  wide  QCL,  was  indeed  115  and 
204  mW  before  and  after  the  deposition  of  the  HR  coating. 
Note  that  the  total  output  power  after  backfacet  coating 
should  be  as  high  as  291  mW  given  the  estimated  collection 
efficiency.  The  above  improvement  was  accompanied  by  a 
decrease  of  approximately  20%  of  the  threshold  current  den¬ 
sity  of  both  devices.  From  this  change,  an  estimate  of  the 
waveguide  losses  was  deduced  and  yields  8.3  and  9  cm-1  for 
the  broad  and  the  narrow  laser,  respectively. 

A  remarkable  feature  is  that  both  lasers  work  in  continu¬ 
ous  mode  at  400  K  with  more  than  10  mW  of  output  power. 

It  exceeds  the  previous  maximum  cw  temperature  operation 
of  QCLs  by  approximately  60  K.11  At  room  temperature,  the 
maximum  collected  output  power  reaches  80  and  204  mW 
with  a  corresponding  slope  efficiency  dP/dl  equal  to  372 
and  473  mW/A  for  the  narrow  and  the  broad  laser,  respec¬ 
tively.  In  a  QCL,  as  the  current  is  increased  above  threshold, 
an  electron  creates  a  number  of  photons  equal  to  the  number 
of  stages  N.  It  is  thus  appropriate  to  define  a  differential 
quantum  efficiency  per  stage  as  rjd=e  ( Nhv)~l  dP/dl  with 
hv  the  photon  energy.  Note  that  P  in  this  equation  is  equal  to 

the  power  collected  from  a  single  facet  of  a  HR  coated  laser. 
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If  the  70%  collection  efficiency  of  our  setup  is  taken  into 
account,  pd  yields,  in  the  case  of  the  HR  coated  7.5  gm  wide 
QCL,  13.1%  and  3.32%  at,  respectively,  300  and  400  K. 
With  these  numbers  and  the  estimate  of  the  internal  loss 
^waveguide’  the  internal  quantum  efficiency  (per  stage)  rjj  can 
be  found  from  the  standard  formula  rjt=  rjd(\ 
+  “waveguide/ “mirror) >  where  amirror=2.15  cm-1  is  the  mirror 
loss.  When  rjd=  13.1%,  the  internal  quantum  efficiency  is 
about  63.7%,  which  is  comparable  to  standard  edge-emitting 
diode  lasers. 

The  value  of  the  threshold  current  density  y'th  measured 
in  cw  and  pulsed  mode  (1.5%  duty  cycle,  80  KHz  repetition 
rate)  is  plotted  as  a  function  of  temperature  in  Fig.  2(a).  The 
pulsed  threshold  current  of  the  narrow  device  is  approxi¬ 
mately  44%  larger  than  the  value  found  for  the  broad  laser. 
This  percentage  difference  is  essentially  constant  over  the 
entire  temperature  range  and  is  explained  by  the  reduction  of 
the  overlap  of  the  optical  mode  with  the  gain  medium  as  the 
stripe  width  narrows.  jth  is  indeed  proportional  to  the  total 
optical  losses  divided  by  the  overlap  factor  T.  Using  the 
waveguide  loss  mentioned  earlier  in  the  text  and  the  values 
of  T  calculated  with  a  two-dimensional  (2D)  beam- 
propagation  modeling  program  (38%  and  55%  for  the  nar¬ 
row  and  the  broad  QCL,  respectively),  jth  is  found  to  be  54% 
larger  for  the  3  gm  wide  lasers  compared  to  the  broad  de¬ 
vice,  in  good  agreement  with  the  experimental  data.  The  cw 
threshold  current  densities  of  the  3  gm  wide  laser  are 

1.9  kA/cm2  at  room  temperature  and  3.5  kA/cm2  at  400  K. 
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FIG.  3.  Room  temperature  cw  differential  resistance  of  the  7.5  jmm  HR 
coated  device  as  a  function  of  current  after  8,  200,  and  500  h  of  cw  opera¬ 
tion  at  370  K.  The  discontinuity  indicates  the  position  of  the  threshold  cur¬ 
rent.  Inset:  wall-plug  efficiency  as  a  function  of  current  for  temperatures 
ranging  from  240  to  400  K. 

As  in  pulsed  mode,  the  corresponding  values  for  the  wider 
laser  are  lower  (1.39  kA/cm2  at  300  K  and  still  less  than 
3  kA/cm2  at  400  K). 

It  is  apparent  from  Figs.  1  and  2(a)  that  the  7.5  /um  wide 
ridge  is  the  most  sensitive  of  the  two  samples  with  respect  to 
temperature  changes.  The  characteristic  temperature  T0 , 
which  describes  phenomenologically  the  increase  of  y'th  with 
temperature,  has  a  similar  value  in  pulsed  mode  for  both  the 
narrow  and  the  broad  laser.  A  much  larger  T0  is  found,  how¬ 
ever,  in  continuous  mode  for  the  3  jmm  wide  device,  which 
confirms  that  a  narrow  stripe  width  leads  to  performance  less 
sensitive  to  temperature.  This  conclusion  is  supported  by  Fig. 
2(b),  in  which  the  slope  efficiency  is  plotted  as  a  function  of 
temperature.  The  value  of  dP / dl  clearly  decreases  faster  for 
the  device  with  a  broader  ridge  as  the  temperature  increases 
and  above  380  K,  the  slope  efficiency  of  the  narrowest 
sample  even  becomes  larger.  The  better  heat  removal  capa¬ 
bility  of  the  narrow  stripe  laser  translates  into  a  higher  value 
of  the  thermal  conductance  Gth=(RthA)_1.  Here  A  is  the  de¬ 
vice  area  and  Rth  (10.84  and  11.79  K/W  for  the  3  and 
7.5  jmm  device  at  300  K)  is  the  thermal  resistance  which  can 
be  determined  from  the  threshold  current  density  versus  tem¬ 
perature  curves  shown  in  Fig.  2(b).6  At  300  K,  the  thermal 
conductance  yields  942  and  394  W/Kcm2  for,  respectively, 
the  narrow  and  the  broad  device.  Note  that  the  Gth  ratio  of 
2.39  is  very  close  to  the  stripe  width  ratio  of  2.5.  The  larger 
value  of  Gth  found  for  the  3  jmm  stripe  is  understandable  as 
the  ratio  of  volume  to  area  is  more  favorable  in  terms  of  heat 
dissipation  for  narrower  ridges.  Similar  conclusions  were 
reached  in  previous  studies.1 

The  inset  of  Fig.  2(b)  shows  typical  emission  spectra 
obtained  above  threshold  with  the  devices  investigated  in  the 
present  letter.  The  different  curves  were  measured  at  300  K 
with  a  3  jmm  wide  laser  operating  in  continuous  mode.  The 
emission  wavelength  is  8.38  jmm  and  little  shorter  than  the 
one  of  Refs.  10.  Note  that  the  optical  spectra  broadened  and 
showed  lasing  in  many  modes  when  the  devices  were 
pumped  at  current  densities  at  least  30%  above  jth.  A  detailed 
study  of  this  phenomenon  will  be  published  elsewhere.  The 
wall-plug  efficiency,  defined  as  the  total  optical  output  power 
divided  by  the  electrical  input  power,  is  shown  in  the  inset  of 
Fig.  3  for  the  HR  coated  7.5  jmm  wide  device.  The  I-V  and 
L-I  curves  presented  in  Fig.  1(b)  were  used  in  the  calcula- 
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tions.  The  wall-plug  efficiency  is  as  high  as  4.6%  at  240  K 
and  0.18%  at  400  K.  Given  the  estimated  collection  effi¬ 
ciency  of  our  setup,  these  values  should  be  approximately  a 
factor  of  1/0.7  larger. 

Preliminary  results  obtained  in  the  context  of  accelerated 
aging  measurements  are  also  promising.  The  7.5  jmm  wide 
laser,  whose  performance  is  described  in  this  letter  was  op¬ 
erated  for  long  periods  of  time  in  cw  at  370  K  and  a  constant 
current  of  630  mA,  corresponding  to  an  output  power  of 
30  mW  and  a  current  density  of  2.8  kA/cm2.  At  regular  in¬ 
tervals,  the  cw  L-I  and  I-V  curves  were  measured  at  300  K. 
After  more  than  500  h,  only  negligible  variations  in  the 
threshold  current  and  slope  efficiency  were  observed.  Figure 
3  shows,  for  example,  the  differential  resistance  of  the  de¬ 
vice,  deduced  from  the  I-V  curves  measured  at  room  tem¬ 
perature  after  8,  200,  and  500  h  of  cw  operation  at  370  K. 
The  threshold  current  corresponds  to  the  discontinuity  in 
each  of  the  traces.  It  varies  between  curves  by  less  than 
2  mA.  This  change  (approximately  0.6%)  is  insignificant  as 
it  is  comparable  to  the  measurement  error.  The  latter  is  equal 
to  3  mA  and  is  due  to  the  fact  that  the  heat  sink  temperature 
is  only  controlled  within  0.2  K.  Note  that  the  box  in  which 
the  laser  was  mounted  was  not  purged  with  dry  air  or  nitro¬ 
gen.  A  more  detailed  description  of  the  results  will  be  given 
elsewhere. 

In  conclusion,  8.38  jmm  wavelength,  HR  coated  QC  la¬ 
sers  grown  by  low-pressure  MOVPE  and  processed  into  nar¬ 
row  buried  heterostructures  were  operated  in  cw  mode  up  to 
400  K.  High  cw  output  powers  as  high  as  204  mW  at  room 
temperature  and  more  than  10  mW  at  400  K  were  measured. 
The  laser  characteristics  reported  here  represent  a  significant 
improvement  in  the  performance  of  MOVPE-grown  QCLs. 
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The  authors  report  the  fabrication  of  high-power  strained  quantum  cascade  lasers  working  in 
continuous  mode  above  370  K.  The  devices,  processed  in  narrow  buried  heterostructures,  were 
grown  by  low-pressure  metal  organic  vapor-phase  epitaxy.  Continuous  wave  output  power  as  high 
as  312  mW  at  300  K  was  obtained  at  a  wavelength  of  5.29  /xm  from  a  3.25  mm  long,  7.5  /xm  wide 
laser  with  a  high-reflectivity  back  facet  coating.  The  slope  efficiency  was  in  excess  of  1.5  W/A  and 
the  power  conversion  efficiency  reached  almost  5%.  ©  2006  American  Institute  of  Physics. 
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In  the  past  few  years,  the  performance  of  quantum  cas¬ 
cade  lasers1  (QCLs)  grown  by  molecular  beam  epitaxy 
(MBE)  has  improved  tremendously  and  high-power  continu¬ 
ous  wave  (cw)  operation  at  room  temperature  has  been  since 
routinely  achieved  at  various  wavelengths.2,3  Short- 
wavelength  (5.25  yam)  QCLs  with  a  maximum  cw  output 
power  of  480  mW  were  reported  in  Ref.  4.  This  result  was 
obtained  at  room  temperature  (298  K)  with  a  strained  QC 
structure  based  on  a  double-phonon  resonance1  processed 
into  11  /xm  wide  ridge  waveguide.  Together  with  the  large 
cw  output  power  achievable  by  QCLs  at  room  temperature, 
the  wavelength  range  accessible  to  these  devices  makes  them 
very  attractive,  in  particular,  for  compact  gas  sensors  based 
on  photoacoustic  spectroscopy.5 

Metal  organic  vapor-phase  epitaxy  (MOVPE)  is  a  widely 
established  platform  for  the  high- volume  production  of  low 
cost  and  reliable  semiconductor  lasers.  This  technique  allows 
the  deposition  of  uniform  and  thick  layers  at  high  rates.  It 
offers  also  excellent  stability  over  long  growth  runs  and  short 
machine  downtime.  MOVPE  is  thus  particularly  well  suited 
for  the  industrial  production  of  QCLs.  Recent  work  showed 
that  the  performance  obtained  with  lattice-matched  and 
strained  QCLs  deposited  by  MOVPE  is  comparable  with  that 
of  similar  structures  grown  by  MBE.6-8  Recently,  we  dem¬ 
onstrated  high-power  (=200  mW)  cw  operation  at  room 
temperature  of  X  =  8.38  jam  lattice-matched  QCLs  grown  by 
MOVPE.9 

In  the  present  letter,  we  report  the  cw  operation  of 
strained  QCLs  at  temperatures  ranging  from  300  to  385  K. 
The  devices  studied  were  processed  into  buried  heterostruc¬ 
tures  and  a  high-reflectivity  (HR)  coating  was  evaporated  on 
the  back  facet.  With  a  7.5  /xm  wide  3.25  mm  long  back  facet 
coated  laser,  we  measured  low  threshold  current  densities  at 
room  temperature  equal  to  1.23  kA/cm2  in  pulsed  mode  and 
to  1.53  kA/cm2  in  cw,  together  with  high  cw  output  power  in 
excess  of  300  mW  at  300  K.  The  maximum  operating  tem- 
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perature  was  as  high  as  376  K  for  this  particular  device.  This 
represents  a  major  improvement  compared  to  our  previous 
results  obtained  with  similar  devices,  which  operated  only  up 
to  approximately  280  K  in  cw  mode.8 

The  QCL  structure  studied  in  the  present  work  was  de¬ 
posited  on  a  heavily  doped  n- type  (xz~  3  X  1018  cm3)  InP 
substrate  by  low-pressure  (76  Torr)  MOVPE  without  growth 
interruption  at  the  barrier/quantum  well  interfaces.  The 
growth  conditions  were  similar  to  the  ones  described  in  Refs. 
7-9.  The  deposition  rate  was  0.5  nm/s  throughout  the  struc¬ 
ture,  except  for  the  active  region,  which  was  grown  at  a 
slower  rate  of  0.1  nm/s.  The  cladding  layers  consist  of  a 
3.5  /xm  thick  InP  layer  doped  with  silicon  (xz~  1 
X  1017  cm3).  The  lower  cladding  layer  helps  in  decoupling 
the  optical  mode  from  the  heavily  doped  substrate, 
preventing  substantial  free-carrier  absorption.  The  wave¬ 
guide  core  is  composed  of  two  0.33  /xm  thick  n-doped  {n 
~  3  X  1016  cm3)  InGaAs  layers  and  the  active  region.  The 
latter  comprises  30  stages  and  is  based  on  a  double-phonon 
resonance  design  identical,  except  for  the  doping  level  in 
the  injector,  to  the  one  reported  by  Hofstetter  et  al .10  The 
growth  ended  with  a  0.5  jam  thick  InP  plasmon-enhanced 
confinement1  layer  doped  with  silicon  (xz~5  X  1018  cm3), 
followed  by  two  heavily  doped  contact  layers  (10  nm, 
n  +  InP/ 20  nm,  n+  InGaAs).  Overall,  the  structure  de¬ 
scribed  above  is  very  similar  to  the  double-phonon  laser  re¬ 
ported  in  a  previous  article.8  The  main  difference  resides  in 
the  doping  level  being  lowered  by  a  factor  of  2  in  ah  of  the 
doped  layers  comprising  the  claddings  and  waveguide  core 
(including  the  active  region). 

The  grown  structure  was  processed  into  buried  hetero¬ 
structure  lasers.  Ridges  3  and  7.5  /xm  wide  were  defined  by 
conventional  photolithography  techniques  and  etched  down 
to  approximately  half  of  the  lower  cladding  layer  by  reactive 
ion  etching.  Using  MOVPE,  InP  doped  with  Fe  was  then 
selectively  regrown  around  the  ridges  to  enhance  the  lateral 
heat  dissipation  from  the  active  region.  The  level  of  Fe  dop¬ 
ing  was  chosen  to  slightly  exceed  the  background  doping  of 
the  MOVPE  reactor.  Since  Fe  atoms  introduced  in  InP  act  as 
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FIG.  1.  (a)  V-I/L-I  curves  obtained  at  different  temperatures  in  cw  mode 
with  a  7.5  jmm  wide,  3.25  mm  long  HR-coated  device,  (b)  Threshold  current 
density  measured  in  pulsed  and  cw  modes  vs  temperature.  The  experimental 
points  were  fitted  with  the  empirical  relation  jth=j0  exp(T/r0). 

deep  carrier  traps,  there  are  essentially  no  free  electrons  in 
the  regrown  InP  layer,  resulting  in  insulating  and  optically 
transparent  regions  around  the  laser  ridges.  The  next  fabrica¬ 
tion  steps  include  the  evaporation  of  Ti/Au  top  contacts  fol¬ 
lowed  by  the  electroplating  of  a  5  jmm  thick  Au  layer.  Thin¬ 
ning  of  the  substrate  and  the  deposition  of  a  Ge/Au  back 
contact  finished  the  processing.  The  devices  were  subse¬ 
quently  cleaved  into  cavities  of  different  lengths  before  being 
soldered  episide  up  onto  Cu  heat  sinks.  A  HR  coating  con¬ 
sisting  of  A1203/Au  (200/30  nm)  layers  was  evaporated  on 
the  back  facet  of  the  QCLs.  The  voltage  and  output  power 
versus  current  (V-I/L-I)  characteristics  were  obtained  with 
the  devices  placed  into  a  box  equipped  with  thermoelectric 
cooler.  The  output  power  was  measured  using  a  polished 
metallic  pipe  guiding  the  light  from  the  laser  front  facet  di¬ 
rectly  onto  the  surface  of  a  calibrated  thermopile  detector. 
The  collection  efficiency  of  this  setup  is  considered  to  be 
near  100%.  The  spectral  characterization  was  performed  with 
a  Fourier  transform  infrared  spectrometer  equipped  with  a 
deuterated  triglycine  sulphate  detector. 

The  V-I/L-I  curves  obtained  in  continuous  mode  from  a 
3.25  mm  long,  7.5  jmm  wide  and  HR-coated  QCLs  are 
shown  in  Fig.  1(a).  The  cw  threshold  current  density  jth  for 
this  device  was  1.53  kA/cm2  at  room  temperature  and  the 
corresponding  slope  efficiency  was  1.570  W/A.  The  largest 
output  powers  obtained  at  300  and  360  K  were  312  and 

57  mW,  respectively.  The  highest  cw  temperature  operation 
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rmax  was  376  K  and  by  reducing  the  temperature  to  260  K, 
more  than  500  mW  of  cw  output  power  was  achieved.  The 
maximum  wall-plug  efficiency,  defined  as  the  total  optical 
power  divided  by  the  electrical  power,  reached  4.6%  at 
300  K.  This  device  has  performance  comparable  to  the  best 
and  most  recent  results  obtained  with  slightly  broader  MBE- 
grown  QCLs  emitting  at  the  same  wavelength.4 

Figure  1(b)  is  a  plot  of  the  threshold  current  densities 
obtained  at  various  heat  sink  temperatures  from  the  HR- 
coated  device  described  in  the  previous  paragraph.  These 
data  were  deduced  from  the  cw  measurements  presented  in 
Fig.  1(a)  and  from  measurements  obtained  in  pulsed  mode. 
In  the  latter  case,  the  repetition  rate  was  80  kHz  and  the 
pulse  width  was  125  ns.  From  these  data,  high  characteristic 
temperatures  T0  equal  to  205.2  K  in  pulsed  mode  and 
169.3  K  in  cw  mode  were  deduced.  T0  describes  empirically 
the  increase  in  threshold  current  versus  temperature  accord¬ 
ing  to  an  exponential  function  jth(r)=70  exp (77 To).  The  ther¬ 
mal  resistance  of  the  device  Rth  can  also  be  calculated  from 
the  measurements  presented  in  Fig.  1(b)  and  yields 
12.14  K/W  at  300  K.11  It  corresponds  to  a  thermal  conduc¬ 
tance  Gth=l/(ARth)  =  337.9  W/(K  cm2),  where  A  is  the  de¬ 
vice  area.  Using  the  value  found  for  7?th,  the  temperature 
difference  A  T  between  the  heat  sink  and  the  active  region 
reaches  at  300  K  AT=RthFelec=42.8  K.  Here  Felec  is  the 
electrical  power  dissipated  in  the  device. 

Figures  2(a)  and  2(b)  show  data  qualitatively  similar  to 
those  displayed  in  Fig.  1.  This  second  set  of  measurements 
was,  however,  obtained  from  a  3  jmm  wide  laser.  The  other 
device  characteristics  such  as  the  cavity  length  were  other¬ 
wise  identical.  The  cw  output  power  was  as  high  as  97  mW 
at  300  K,  which  is  noteworthy  given  the  narrow  width  of  this 
QCL.  The  corresponding  slope  efficiency  was  1.19  W/A. 
The  threshold  current  densities  measured  at  room  tempera¬ 
ture  were  1.39  and  1.43  kA/cm2  in  pulsed  and  cw  modes, 
respectively.  The  value  of  jth  in  pulsed  mode  is  only  slightly 
larger  than  the  one  obtained  with  the  7.5  jmm  wide  structure. 
This  is  consistent  with  the  results  of  two-dimensional  mode 
profile  simulations,  which  showed  that  the  mode  overlap  fac¬ 
tor  T  is  comparable  for  both  the  narrow  and  broad  samples 
(53.3%  and  61.6%,  respectively).  The  inset  of  Fig.  3  shows  a 
typical  spectrum  obtained  at  a  current  5%  above  threshold 
with  the  narrow  QCL.  Note  that  when  the  devices  studied  are 
pumped  30%  or  more  above  jth,  the  optical  spectrum  broad¬ 
ens  and  is  composed  by  many  neighboring  longitudinal 
modes.  A  detailed  study  of  this  phenomenon  will  be  given 
elsewhere.  The  emission  wavelength  is  5.29  jmm  and  is  very 
close  to  the  designed  value.10 

The  thermal  resistance  for  the  narrowest  QCL  is 
7.01  K/W,  which  translates  into  Gth=1461  W/Kcm2.  This 
value  is  comparable  to  the  one  found  for  state-of-the-art  bur¬ 
ied  heterostructure  QCLs  whose  performance  was  reported 
in  Ref.  12.  The  ratio  of  the  thermal  conductance  found  for 
the  broad  and  narrow  samples  does  not  match  the  ratio  of  the 
lasers  widths.  However,  it  indicates  clearly  that  lateral  heat 
transfer  is  a  major  factor  in  the  thermal  resistance  of  the 
device.  This  is  also  in  agreement  with  the  fact  that  the  quan¬ 
tity  T0  has  a  comparable  value  in  pulsed  and  cw  modes  for 
the  narrow  device  (pulsed  mode,  T0=  189.3  K;  cw  mode, 
T0=  185.7  K),  while  T0  in  the  case  of  the  7.5  jmm  QCL  de¬ 
creases  by  almost  35  K  (pulsed  mode,  T0  =  205.2  K;  cw 
mode,  T0=  169.3  K). 
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FIG.  2.  (a)  V-I/L-I  curves  obtained  at  different  temperatures  in  cw  mode 
with  a  3  /mm  wide,  3.25  mm  long  HR-coated  device,  (b)  Threshold  current 
density  measured  in  pulsed  and  cw  modes  vs  temperature.  The  experimental 
points  were  fitted  with  the  empirical  relation  jth=j0  exp(T/T0). 

Three  jum  wide,  HR-coated  devices  with  different  cavity 
lengths  were  also  measured.  The  longest  sample  (4.25  mm) 
had  the  highest  cw  operating  temperature  (385  K)  and  the 
largest  output  power,  which  reached  136  mW  at  300  K.  The 
threshold  current  density  measured  in  pulsed  mode  for  1.09, 
1.75,  2.25,  3.25,  and  4.25  mm  long  cavity  is  reported  in  Fig. 
3  as  a  function  of  the  inverse  of  the  cavity  length.  From  these 
data,  the  waveguide  losses  aw  were  deduced  assuming  that 
the  gain  cross  section  does  not  depend  on  the  current 
density.13  The  value  of  aw  calculated  from  the  fit  of  the  ex¬ 
perimental  points  is  only  3.8  cm-1,  which  compares  well 
with  the  numbers  reported  in  Ref.14 

In  summary,  we  have  presented  the  high-power  cw  op¬ 
eration  at  high  temperature  of  5.29  jmm  wavelength  strained 
QCLs  grown  by  low-pressure  MOVPE.  The  characteristics 
of  our  lasers  are  comparable  to  the  best  and  most  recent 
results  obtained  from  MB E- grown  material.  Our  devices 
have  the  highest  operating  cw  temperature  reported  so  far  for 
QCLs  at  this  wavelength.  The  present  work  demonstrates 
that  MOVPE  meets  all  the  requirements  necessary  to  grow 


FIG.  3.  Pulsed  threshold  current  density  and  device  thermal  resistance  mea¬ 
sured  at  300  K  with  3  jmm  wide  HR-coated  lasers  having  a  different  cavity 
length.  The  inset  shows  an  optical  spectrum  measured  with  the  3  gm  wide, 
3.25  mm  long  HR-coated  device  operated  in  cw  mode  close  to  threshold. 

high-performance  state-of-the-art  QCLs,  even  with  demand¬ 
ing  material  constraints  such  as  large  strain. 
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We  report  the  observation  of  a  coherent  multimode  instability  in  quantum  cascade  lasers  (QCLs),  which  is 
driven  by  the  same  fundamental  mechanism  of  Rabi  oscillations  as  the  elusive  Risken-Nummedal-Graham- 
Haken  (RNGH)  instability  predicted  40  years  ago  for  ring  lasers.  The  threshold  of  the  observed  instability  is 
significantly  lower  than  in  the  original  RNGH  instability,  which  we  attribute  to  saturable-absorption  nonlin¬ 
earity  in  the  laser.  Coherent  effects,  which  cannot  be  reproduced  by  standard  laser  rate  equations,  can  play 
therefore  a  key  role  in  the  multimode  dynamics  of  QCLs,  and  in  lasers  with  fast  gain  recovery  in  general. 
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The  fundamental  coherent  mechanism  that  can  destabilize 
a  single-mode  laser  was  predicted  in  the  early  1960s  [1]  and 
was  later  extended  to  multimode  lasers  [2,3]  where  it  became 
known  as  the  Risken-Nummedal-Graham-Haken  (RNGH) 
instability.  These  instabilities  became  classic  landmarks  for 
the  general  field  of  nonlinear  dynamics  [4,5]  because  they 
emerge  in  conceptually  the  simplest  laser  model,  which  in 
the  single-mode  case  was  shown  to  be  equivalent  to  the  Lor- 
entz  model  of  deterministic  chaos  [6].  Another  feature  that 
makes  these  instabilities  so  interesting  and  unique  is  their 
coherent  nature  that  involves  the  polarization  of  the  medium 
as  a  nontrivial  dynamical  variable.  Most  other  physical 
mechanisms  that  can  drive  a  laser  from  a  single-mode  to  a 
multimode  regime,  such  as  spatial  and  spectral  hole  burning, 
Q  switching,  and  saturable  absorption  [7,8],  can  be  ad¬ 
equately  described  within  the  standard  rate  equation  formal¬ 
ism,  in  which  the  polarization  of  the  active  medium  is  adia- 
batically  eliminated.  Both  the  single  mode  [1]  and  the 
multimode  [2,3]  instabilities  cannot  be  explained  by  the  rate 
equations.  Such  coherent  effects  can  be  only  observed  when 
the  polarization  is  driven  faster  than  or  on  a  time  scale  com¬ 
parable  to  the  dephasing  time  T2  [9]. 

The  origin  of  the  two  coherent  phenomena  mentioned 
above  is  the  oscillation  of  the  population  inversion  at  the 
Rabi  frequency  QRabi  that  takes  place  when  the  intracavity 
laser  intensity  becomes  large.  This  results  in  a  modification 
of  the  gain  spectrum  and  the  emergence  of  sidebands  sepa¬ 
rated  from  the  maximum  of  the  gain  curve  by  an  amount 
corresponding  to  the  Rabi  frequency.  These  sidebands  can  be 
regarded  as  a  manifestation  of  parametric  gain.  The  instabil¬ 
ity  sets  in  when  the  intracavity  power  is  sufficiently  large: 
the  Rabi  angular  frequency  HRabi  has  to  be  greater  than  the 
relaxation  time  scales  of  the  gain  medium  [more  precisely, 
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flRabi  is  sufficiently  greater  than  (r1r2)_1/2,  where  Tx  is  the 
gain  relaxation  time].  The  instability  threshold  is  often  called 
the  second  laser  threshold  due  to  its  universal  nature. 

Pioneering  theoretical  works  stimulated  extensive  experi¬ 
mental  studies  that  finally  resulted  in  the  observation  of  the 
Lorenz-type  chaos  in  a  far-infrared  single-mode  laser  [10]. 
However,  despite  almost  40  years  of  efforts,  the  experimen¬ 
tal  demonstration  of  the  multimode  RNGH  instability  has 
remained  somewhat  controversial  [11-16]. 

In  lasers  with  long  gain  recovery  compared  to  the  cavity 
round-trip  time,  the  instability  caused  by  a  saturable  absorber 
can  often  lead  to  mode  locking  [8].  When  the  gain  recovery 
time  is  short  compared  with  the  cavity  round-trip  time,  it  is 
usually  assumed  that  the  laser  dynamics  becomes  very  primi¬ 
tive  and  uninteresting  (so-called  class  A  laser).  In  this  case 
mode  locking  is  impossible  according  to  conventional  theory, 
and  the  relaxation  oscillation  frequency  becomes  purely 
imaginary  [17].  Surprisingly,  as  we  show  in  this  Rapid  Com¬ 
munication,  it  is  under  these  conditions  that  the  elusive 
RNGH  instability  can  be  observed.  We  show  that  quantum 
cascade  lasers  (QCLs)  are  uniquely  suited  for  studying  these 
coherent  effects  which,  along  with  spatial  hole  burning 
(SHB),  become  a  key  factor  in  dictating  the  dynamics  of  the 
laser. 

QCLs,  because  they  are  based  on  ultrafast  tunneling  and 
phonon-limited  intersubband  transitions,  belong  to  the  class 
of  lasers  which  have  an  extremely  fast  gain  recovery,  in  the 
range  of  a  few  ps  [18].  Recent  experiments  showed  indeed 
that  the  gain  recovers  within  a  few  ps,  which  is  approxi¬ 
mately  an  order  of  magnitude  shorter  than  the  cavity  round- 
trip  time  [19].  Since  its  invention  in  1994,  QCLs  have  un¬ 
dergone  tremendous  improvement  [20].  Recent  development 
of  low  loss,  high  power  QCLs  [21,22]  enables  the  study  of 
previously  under  investigated  aspects,  such  as  the  richness  of 
the  optical  spectrum  and  the  ultrafast  dynamics  of  these  de¬ 
vices.  In  Ref.  [23],  strong  evidence  of  self-pulsations  at  the 
cavity  round-trip  frequency  was  reported  in  QCLs,  in  par- 
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ticular  a  large  broadening  of  the  spectrum  above  the  thresh¬ 
old  of  this  instability  was  observed.  However,  no  detailed 
pulse  characterization  was  provided.  The  technological  po¬ 
tential  of  QCLs  calls  for  a  better  understanding  of  the  inter¬ 
play  of  various  instabilities  in  the  parameter  regime  dictated 
by  these  lasers.  Moreover,  the  Rabi  frequency  in  QCLs  at  the 
power  levels  of  a  few  hundred  mW  is  of  the  order  of  a  few 
THz,  much  larger  than  the  spacing  of  Fabry-Perot  modes. 
Therefore  coherent  effects  should  be  easily  observable  in 
QCLs. 

In  this  Rapid  Communication  we  present  a  clear  experi¬ 
mental  demonstration  of  a  coherent  instability,  driven  by  the 
same  mechanism  as  the  RNGH  instability.  It  is  identified  in 
the  most  direct  manner,  by  demonstrating  in  the  optical  spec¬ 
trum  of  QCLs  a  splitting  corresponding  to  twice  the  Rabi 
frequency. 

The  instability  observed  differs  in  some  aspects  from  the 
original  RNGH  instability  [2,3].  The  threshold  of  instability 
can  be  as  low  as  a  few  tens  of  percent  above  the  laser  thresh¬ 
old,  as  shown  in  Fig.  1(a).  In  addition,  the  pure  RNGH  in¬ 
stability  typically  gives  rise  to  spectra  with  one  central  mode 
and  two  sidebands  separated  from  it  by  the  Rabi  frequency, 
whereas  in  our  experiments  we  observed  two  peaks  only,  x 

similarly  to  Ref.  [11].  However,  the  mechanism  of  the  insta-  “ 

bility  is  the  same  in  essence,  namely  the  Rabi  oscillations  of  c 
the  population  inversion  due  to  coherent  effects.  The  differ¬ 
ences  from  the  RNGH  instability  as  it  occurs  in  ideal  condi¬ 
tions  [2,3]  can  be  attributed  to  the  presence  of  saturable  ab¬ 
sorption  and  SHB. 

The  QCLs  studied  were  fabricated  from  two  different  wa¬ 
fers  (wafer  nos.  3251  and  3252)  grown  by  metalorganic  va¬ 
por  phase  epitaxy.  The  devices  were  processed  into  buried 
hetero structure  lasers,  in  which  an  insulating  Fe-doped  InP 
layer  is  regrown  after  etching  of  the  ridges  [21,22].  The  ac¬ 
tive  region  of  all  the  samples  tested  is  based  on  a  four- 
quantum- well  design,  which  relies  on  a  double  phonon  reso¬ 
nance  to  achieve  population  inversion  [24].  Note,  however, 
that  the  multimode  operation  described  in  the  present  paper 
was  also  observed  with  lasers  based  on  so-called  three- 
quantum- well  designs  [18].  Figure  1(a)  shows  the  optical 
spectra  of  a  laser  operated  in  continuous  wave  (cw)  at  room 
temperature.  The  active  region  of  this  laser  is  3  jum  wide  and 
its  emission  wavelength  is  close  to  8.47  gm  (wafer  no. 

3251).  The  laser  was  cleaved  into  a  3-mm-long  bar  and  sol¬ 
dered  with  indium  onto  a  copper  heat  sink.  The  spectra  were 
measured  by  a  Fourier  transform  infrared  spectrometer 
(FTIR)  equipped  with  a  deuterated  triglycine  sulphate 
(DTGS)  detector. 

As  shown  in  Fig.  1(a),  the  laser  spectrum  is  single  mode 
close  to  threshold  and  broadens  as  the  pumping  current  in¬ 
creases,  splitting  into  two  separated  humps.  The  difference 
between  the  weighted  centers  of  the  two  peaks  increases  lin¬ 
early  as  a  function  of  the  square  root  of  the  collected  output 
power  from  one  facet,  as  shown  in  Fig.  1(b)  (square  dots 
with  the  dashed  line  as  its  best  fit).  The  Rabi  angular  fre¬ 
quency  flRabi  can  be  easily  calculated  using  the  formula 
HRabi = ijlE Ih=/JL^ Inl^J  ( ce)lfi ,  where  /z  is  the  electron 
charge  times  the  matrix  element  of  the  laser  transition 
(=2.54  nm).  7ave  is  the  average  intracavity  intensity  in  the 
gain  region,  which  can  be  derived  from  the  measured  output 
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FIG.  1.  (Color  online)  (a)  Optical  spectra  vs  pumping  ratio 
0  / 7th)  above  threshold  obtained  in  cw  at  300  K  with  a  3-/zm-wide 
buried  heterostructure  lasers  emitting  at  8.47  /zm.  For  1  <///th 
<1.2  the  spectra  are  identical  to  j/j±=  1.2.  (b)  Spectral  splitting 
and  twice  the  Rabi  frequency  HRabi/(27r)  vs  square  root  of  output 
power  collected  from  a  single  laser  facet.  The  different  quantities 
reported  on  the  graph  were  deduced  from  the  experimental  data 
shown  in  (a).  The  dashed  line  is  a  least- square  linear  fit  of  the  data. 

power,  c  is  the  speed  of  light  in  vacuum  and  n  is  the  back¬ 
ground  refractive  index.  For  all  the  values  of  the  intensity 
corresponding  to  the  spectra  reported  in  Fig.  1(a),  nRabi  was 
calculated,  multiplied  by  a  factor  2  and  then  added  to  Fig. 
2(b)  (solid  line).  A  very  good  agreement  is  found  between 
the  experimental  splitting  and  twice  the  estimated  Rabi  fre¬ 
quency.  Both  curves  fall  indeed  well  within  the  error  bars 
[25].  As  mentioned  before,  the  theory  behind  the  RNGH  in¬ 
stability  predicts  that  the  large  intracavity  intensity  will  result 
in  parametric  gain  at  frequencies  detuned  from  the  maximum 
of  the  gain  curve  by  the  Rabi  frequency.  The  agreement  men¬ 
tioned  above  is  thus  a  strong  indication  of  the  RNGH  insta¬ 
bility  in  QCLs. 

In  order  to  better  understand  the  experimental  spectra  of 
the  QCLs  presented  in  Fig.  1(a),  we  use  a  simple  model 
based  on  the  standard  one-dimensional  Maxwell-Bloch  equa¬ 
tions  [9],  where  the  active  medium  is  described  by  an  “open” 
two-level  system  [26].  However,  contrary  to  the  standard 
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FIG.  2.  Results  of  numerical  simulations  of  the  spectra  based  on 
the  Maxwell-Bloch  equations  including  a  saturable  absorber  and 
spatial  hole  burning  for  different  values  of  the  current  density  nor¬ 
malized  to  the  threshold  value. 

unidirectional  Maxwell-Bloch  equations,  we  allow  the  elec¬ 
tromagnetic  field  to  propagate  in  both  directions.  The  waves 
traveling  in  the  two  directions  are  coupled,  as  they  share  the 
same  gain  medium.  This  gives  rise  to  SHB  [7]:  The  standing 
wave  formed  by  a  cavity  mode  imprints  a  grating  in  the  gain 
medium  through  gain  saturation.  As  a  result,  other  modes 
become  more  favorable  for  lasing,  and  a  multimode  opera¬ 
tion  is  triggered. 

In  the  slowly  varying  envelope  approximation,  the  equa¬ 
tions  read 

n  kNfiT  1 

-d,E±  =  T  dzE±  -  i—^j  V±  -  -t(E+,E_)E±,  (la) 

c  2e0n  2 

d,V±  =  “t(A0£±  +  (lb) 

2  n  T2 

#Ao  =  A/2  A°  +  T'(£+77+  +  E-V-  ~  c.c.),  (lc) 
Tx  n 

dAi  =  ±  iT(Kv-  -  V+E-)  ~  •  (Id) 

n  T] 

The  +  and  —  subscripts  label  the  two  directions  of  propaga¬ 
tion.  E  and  77  are  the  slowly  varying  envelopes  of  the  electric 
field  and  the  polarization,  respectively.  The  actual  electric 
field  and  polarization  are  obtained  by  multiplying  E  and  77  by 
el(x)t  (a)  is  the  optical  resonance  frequency)  and  taking  the  real 
part.  The  position-dependent  inversion  is  written  as  the  sum 
of  the  three  terms,  A0,  A2e2lkz,  and  \2e~llkz,  where  (A£)* 
=  A2.  The  inversion  is  thereby  represented  by  two  slowly 
varying  functions  (A0  and  A2),  and  e±2lkz  gives  the  fast  varia¬ 
tion  in  space.  All  the  quantities  mentioned  so  far  are  func¬ 
tions  of  space  z  and  time  t. 
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€  (E+ ,  E_ )  is  the  loss  in  the  cavity  (not  including  the  mirror 
loss),  which  is  allowed  to  be  nonlinear  and  dependent  on  the 
intensity.  In  this  work  we  assume 

€(£+,£_)  =  €0-y(|£+|2+|£_|2),  (2) 

where  €0  is  the  linear  loss  and  y  is  the  self-amplitude  modu¬ 
lation  coefficient  characterizing  the  nonlinear  (saturable)  part 
of  the  loss.  Such  a  saturable  absorption  mechanism  can  come 
from  Kerr  lensing  [8,23],  caused  by  a  nonlinear  refractive 
index  n2  in  the  active  region.  As  the  intensity  increases,  the 
mode  is  more  confined  in  the  plane  transverse  to  the  propa¬ 
gation  direction,  and  the  net  gain  it  experiences  is  greater. 
The  reason  is  twofold:  First,  the  mode  overlaps  more  with 
the  active  region,  leading  to  a  larger  modal  gain  (this  mecha¬ 
nism  is  often  called  “soft  Kerr  lensing”  [27]).  Second,  the 
overlap  with  the  metal  contacts  is  smaller,  leading  to  smaller 
losses  [23]. 

E+  and  E_  satisfy  the  boundary  conditions  E+=rE_  at  the 
z  =  0  boundary  and  rE+=E_  at  the  z=L  boundary  (L  is  the 
cavity  length  and  r~0.53  is  the  reflection  coefficient).  The 
other  quantities  in  Eq.  (1)  are  constants:  k ,  N ,  and  T  are  the 
wave  number  (in  the  material)  associated  with  the  resonance 
optical  frequency,  the  electron  density  in  the  active  region, 
and  the  overlap  factor  between  the  laser  mode  and  the  active 
region,  respectively. 

Figure  2  shows  spectra  that  were  obtained  by  solving  nu¬ 
merically  Eqs.  (1)  with  the  following  parameters:  for  the 
saturable  absorber,  we  used  y=10-8  p,  obtained  from  two- 
dimensional  mode  simulations,  assuming  a  772~10_9^y- 
[23].  The  index  change  due  to  this  n2  at  typical  intracavity 
intensities  is  about  10-3.  The  other  parameters  are  €0 
=  5  cm-1,  7\  =  0.5  ps  [19],  T2=0.07  ps  (corresponding  to  a 
gain  full  width  at  half  maximum  bandwidth  of  4.8  THz),  L 
=  0.3  cm,  and  n-  3,  which  are  typical  values  for  these  lasers. 
N  and  T  are  not  needed  as  long  as  the  pumping  is  expressed 
relative  to  the  lasing  threshold.  Note  that  the  simulated  spec¬ 
tra  presented  in  Fig.  2  are  averaged  over  about  a  microsec¬ 
ond.  Only  then  does  the  average  spectrum  reach  a  steady 
state  and  a  clear  pattern  shows  up.  The  averaging  is  moti¬ 
vated  by  the  fact  that  experimentally  the  spectra  are  acquired 
over  an  even  much  longer  time  scale.  The  envelopes  of  the 
spectra  show  two  clear  peaks  whose  separation  compares 
well  with  twice  the  Rabi  frequency,  similarly  to  the  experi¬ 
ment. 

The  lowering  of  the  RNGH  instability  threshold  by  a  satu¬ 
rable  absorber  can  be  established  analytically  by  means  of 
linear  stability  analysis.  We  propose  this  mechanism  as  the 
main  reason  for  the  observation  of  the  RNGH  instability  at 
much  lower  pumping  than  RNGH  theory  predicts.  In  order  to 
support  this  idea,  we  now  present  spectra  from  another  de¬ 
vice  similar  to  the  one  described  previously.  The  only  differ¬ 
ence  between  the  two  lasers  is  a  shorter  optical  wavelength 
(5.25  /mi)  (wafer  no.  3252)  and  a  wider  active  region 
(5  /mi).  The  two-dimensional  waveguide  simulations  indi¬ 
cate  a  much  weaker  Kerr-lensing  effect  in  these  QCLs  (y  is 
smaller  by  a  factor  of  4),  due  to  the  much  larger  ratio  of 
active  region  width  to  wavelength.  The  measured  optical 
spectra  obtained  at  300  K  in  cw  mode  with  the  5 -/mi  device 
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FIG.  3.  Optical  spectra  vs  pumping  ratio  above  threshold  ob¬ 
tained  in  cw  at  300  K  with  a  5-yum-wide  buried  heterostructure 
lasers  emitting  at  5.25  jam. 

are  shown  in  Fig.  3.  The  data  clearly  show  that  the  laser  is  at 
first  single  mode  close  to  threshold  and  becomes  multimode 
immediately  after  a  slight  increase  of  the  pumping  current. 
The  envelopes  of  the  spectra  consist  of  multiple  peaks,  with 
an  average  separation  0.2  THz,  independently  of  the  pump¬ 
ing.  Numerical  integration  of  Eq.  (1)  without  a  saturable  ab¬ 
sorber  (y=0)  leads  to  spectra  that  qualitatively  agree  with 
the  ones  in  Fig.  3. 

Reference  [15]  suggested  that  the  suppression  of  the  cen¬ 
tral  peak  in  RNGH-type  spectra  can  be  due  to  the  complex 
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level  structure  of  the  gain  medium,  a  dye  molecule  in  that 
case.  We  show  that  SHB  can  also  result  in  the  suppression  of 
the  central  peak  (Fig.  2). 

Our  postulation  of  saturable  absorption  due  to  Kerr- 
lensing  is  supported  by  more  extensive  study  of  different 
devices  beyond  those  shown  in  this  Rapid  Communication. 
First,  we  observed  that  for  the  same  emission  wavelength,  a 
broad  active  region  leads  to  a  less  pronounced  RNGH-type 
signature.  Second,  we  have  also  tested  several  standard  ridge 
waveguide  QCLs,  for  which  the  sidewalls  of  the  ridges  are 
covered  by  the  gold  contact.  For  these  devices  the  coupling 
between  the  optical  mode  and  the  metal  is  expected  to  be 
stronger  and  so  is  the  effect  of  saturable  absorber  due  to  Kerr 
lensing.  The  spectral  behavior  observed  in  this  class  of  de¬ 
vices  is  dominated  by  RNGH-type  instability. 

In  summary,  a  coherent  multimode  instability  in  quantum 
cascade  lasers  (QCLs)  has  been  observed.  It  is  similar 
in  many  ways  to  the  Risken-Nummedal-Graham-Haken 
(RNGH)  instability.  The  threshold  of  the  observed  phenom¬ 
enon  is  significantly  lower  than  in  the  original  RNGH  insta¬ 
bility,  which  is  attributed  to  the  presence  of  a  saturable  ab¬ 
sorber  in  the  laser.  For  devices  with  a  weaker  saturable 
absorber,  the  envelope  of  the  optical  spectrum  consists  of 
many  maxima  whose  separations  are  independent  of  the  int¬ 
racavity  power.  The  nontrivial  shape  of  the  spectrum  can  be 
explained  by  SHB. 
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